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PROJECT  OVERVIEW 


OBJECTIVES 

The  general  objective  of  the  project  is  fundamental  mathematical  modeling  of  a  complex  TSA 
system  with  electrothermal  desorption  step,  with  adsorbers  assembled  of  one  or  more 
activated  carbon  fiber  clot  (ACFC)  cartridge-type,  radial  flow  fixed  beds  and  with  in- vessel 
condensation  of  the  desorbed  vapor.  During  the  first  stage  of  the  project,  a  mathematical 
model  of  a  single  cartridge  without  condensation  was  developed  and  used  for  simulation. 

The  objective  of  this,  second  stage  of  the  project,  is  to  model  the  complete  TSA  system,  for 
adsorbers  with  one  and  more  cartridges.  Besides  the  material  and  energy  balances,  these 
models  have  to  include  the  momentum  balances,  in  order  to  define  the  complex  gas  flow  in 
the  complex  geometry  of  the  adsorber.  The  models  have  to  describe  adsorption, 
electrothermal  desorption  and  electrothermal  desorption  accompanied  with  condensation  of 
the  desorbed  vapor,  as  well  as  the  complete  TSA  cycle.  These  models  will  be  used  for 
prediction  of  the  velocity,  pressure,  concentration  and  temperature  profiles  in  the  adsorbers 
and  calculation  of  the  mass  of  condensed  liquid  and  used  electric  energy.  Applicability  of  the 
developed  mathematical  models  for  optimization  of  the  TSA  process,  with  regard  to  the 
separation  and  purification  factor,  regeneration  of  the  condensed  vapor  and  used  energy,  will 
also  be  investigated. 


STATUS  OF  EFFORT 

Models  for  two  types  of  adsorbers,  one  with  only  one,  and  the  other  with  two  cartridges,  have 
been  developed  in  this  stage.  These  adsorber  types  will  be  called  1 -cartridge  and  2-cartridges 
adsorber,  respectively,  through  this  manuscript.  For  each  adsorber  type,  three  models  were 
built,  in  order  to  describe  three  stages  of  a  compete  TSA  cycle:  adsorption,  electrothermal 
desorption  and  electrothermal  desorption  with  in-vessel  condensation.  These  models  were 
built  using  Femlab,  a  specialized  software  tool  for  modeling  of  complex  systems  with 
complex  geometry.  In  order  to  describe  the  complete  TSA  cycle,  the  models  for  the  three 
stages  were  integrated,  by  using  a  combination  of  Femlab  and  Matlab.  The  models  were 
successfully  used  for  simulation  of  separate  stages  of  the  process  and  of  the  complete  TSA 
cycles,  as  well  as  for  their  optimization. 


ACCOMPLISHMENTS  /NEW  FINDINGS 

The  accomplishments  of  this  project  can  be  perceived  from  two  aspects: 

1 .  The  developed  models  represent  a  mathematical  tool  which  can  be  used  for  simulation 
of  the  investigated  TSA  system,  its  analysis,  predicting  of  the  concentration  and 
temperature  profdes  in  the  system  and  its  optimization.  This  was  the  main  objective  of 
this  project. 

2.  By  developing  Femlab  and  combined  Femlab/Matlab  models  of  this,  very  complex 
system,  involving  a  number  of  coupled  processes  and  phenomena  taking  place  in  a 
complex  geometry,  a  progress  in  the  Femlab  and  Femlab/Matlab  modeling  has  been 
achieved. 
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DETAILED  REPORT 


1.  Introduction 

The  idea  about  regeneration  of  used  adsorbents  by  direct  heating  of  the  adsorbent  particles  by  passing 
electric  current  through  them  (Joule  effect),  was  first  published  in  the  1970s  [1].  Desorption  process 
based  on  this  principle  was  later  named  electrothermal  desorption  [2].  It  was  recognized  as  a 
prospective  way  to  perform  desorption  steps  of  TSA  cycles.  At  the  same  time,  fibrous  activated  carbon 
was  recognized  as  a  very  convenient  adsorbent  form  for  its  realization.  Electrothermal  desorption  has 
some  advantages  over  conventional  methods,  regarding  adsorption  kinetics  and  dynamics  [3,  4]  and 
energy  efficiency  [2,  5].  Some  industrial  applications  of  electrothermal  desorption  have  been  reported 
recently  [6,  7].  Nevertheless,  the  development  of  mathematical  description  of  processes  based  on 
electrothermal  desorption  hasn’t  been  following  the  development  of  the  process  itself,  so  far. 

A  new  TSA  process  with  electrothermal  desorption  step,  based  on  adsorbers  assembled  of  one  or  more 
annular,  cartridge-type,  fixed-beds,  with  in-vessel  condensation,  has  been  presented  recently  [5,  8,  9]. 
Its  schematic  representation  is  given  in  Figure  1. 
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Figure  1.  Overall  sehematie  representation  of  the  ACFC  adsorption  -  rapid  eleetrothermal  desorption  system  (from  Ref  [5]) 


The  aim  of  this  project  is  to  develop  a  rigorous,  fundamental  mathematical  model  of  this  system, 
which  could  be  used  for  simulation,  analysis  and  optimization  of  the  TSA  process  developed  by 
Sullivan  [5]. 


1.1.  Description  of  the  system 

The  adsorbers  of  the  TSA  system  shown  in  Figure  1  are 
composed  of  one  or  more  annular,  radial-flow  cartridge-type 
adsorbent  beds.  A  single  cartridge  is  schematically  shown  in 
Figure  2.  It  is  formed  as  a  cylindrical  roll  of  activated  carbon 
fiber  cloth  (ACFC),  spirally  coiled  around  a  porous  central 
pipe.  The  gas  flow  through  the  adsorber  is  in  the  radial 
direction.  During  the  desorption  step,  electric  current  is 
passed  through  the  activated  carbon  cloth  in  the  axial 
direction,  causing  heat  generation,  heating  of  the  adsorbent 
and  desorption. 
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Figure  2.  Schematic  representation  of  the 
annular  -  radial  flow  -  ACFC  adsorbent  bed 
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The  cartridges  are  situated  in  a  cylindrical  adsorber  vessel.  During  desorption,  the  hot  gas  steam  rich 
with  the  desorbed  vapor  comes  into  contact  with  the  cold  vessel  wall  and  the  vapor  condenses.  In  that 
way,  the  adsorber  vessel  serves  as  a  passive  condenser.  The  bottom  outlet  of  the  adsorber  vessel  is 
funnel-shaped,  in  order  to  facilitate  draining  of  the  condensed  liquid. 
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In  this  phase  of  the  project,  adsorbers  with  one  and  two  cartridges 
will  be  modeled.  The  1 -cartridge  adsorber  is  shown  in  Figure  3  and 
the  2-cartridges  adsorber  in  Figure  4. 

In  Figures  3  and  4  the  directions  of  the  fluid 
flow  are  also  shown,  in  order  to  clarify  the 
processes: 

-  For  the  1 -cartridge  adsorber  the  flow 
directions  during  adsorption  and  desorption 
are  the  same:  the  gas  enters  the  adsorber 
through  the  central  tube,  from  the  bottom,  and, 
after  passing  through  the  adsorbent  bed  in  the 
radial  direction,  exits  from  the  annular  space 
around  the  adsorbent  bed  through  the  funnel- 
shaped  bottom  outlet.  During  desorption,  the 
~  desorbed  vapor  is  partially  condensed  at  the 
cold  outer  wall  of  the  adsorber,  and  the  outlet 
stream  contains  both  gas  and  the  condensed 
liquid. 
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Figure  3.  1 -cartridge 
adsorber 


Figure  4.  2-cartridge  adsorber: 
(A)  adsorption,  (B)  desorption 


-  For  the  2-cartridges  adsorber,  the 
flow  directions  during  adsorption  and 
desorption  are  somewhat  different.  During 
adsorption  the  feed  stream  enters  the  adsorber  from  the  bottom,  through  the 
central  tube  of  the  bottom  cartridge,  and  the  purified  stream  exits  the  adsorber  at  the  top,  through  the 
central  tube  of  the  top  cartridge  (Fig.  4A).  Between  the  inlet  and  the  outlet,  the  stream  flows  through 
both  adsorbent  beds  in  the  radial  direction,  but  in  the  opposite  directions  (through  the  bottom  cartridge 
from  the  center  to  the  periphery,  and  through  the  top  cartridge  from  the  periphery  to  the  center). 
During  desorption,  two  inert  streams  are  introduced  into  the  adsorber,  from  the  bottom  and  the  top, 
through  the  central  tubes  of  both  cartridges,  and  the  gas  and  the  condensed  liquid  exit  the  adsorber 
through  the  funnel  shaped  bottom  outlet  of  the  annular  tube  around  the  cartridges  (Fig  4B).  It  should 
be  noticed  that  the  gas  cannot  flow  directly  between  the  central  tubes  of  the  top  and  the  bottom 
cartridge.  Also,  the  electric  current  supply  for  the  two  cartridges  is  separate. 


The  inlet  gas  flow-rate  during  adsorption  is  generally  several  times  higher  than  during  desorption.  The 
inlet  gas  during  desorption  is  pure  inert. 


1.2.  Brief  survey  of  the  first-year  results 

During  the  first  year,  the  modeling  efforts  were  concentrated  on  modeling  of  a  single  annular, 
cartridge-type,  fixed-bed  adsorber,  without  condensation.  Two  deterministic  models  were  postulated: 
the  first  one  assuming  uniform  adsorbent  density  throughout  the  adsorbent  bed,  and  the  second  one 
taking  into  account  the  structure  of  the  bed  made  of  layers  of  ACFC.  Both  models  were  obtained  as 
complex  sets  of  coupled  nonlinear  PDFs,  DDEs  and  algebraic  equations.  Numerical  procedures  for 
solving  the  model  equations  were  established  and  used  for  simulation  of  adsorption,  electroresistive 
heating,  electrothermal  desorption  and  consecutive  adsorption-desorption,  using  Matlab  software. 
Similar  results  were  obtained  with  the  two  models.  The  results  of  this  stage  are  presented  in  the  Final 
Performance  Report  for  project  FA8655-03-1-3010  [10]. 


2 


1.3.  Outline  of  the  second  year  project 

The  TSA  system  defined  in  Ref.  [5]  is  very  complex,  regarding  both  its  geometry  and  the  processes 
involved.  For  that  reason,  a  decision  was  made  to  perform  the  modeling  of  the  complete  system  by 
using  some  specialized  modeling  software.  Our  choice  was  to  use  FEMLAB,  relatively  new  modeling 
software  based  on  the  finite-element  method  for  solving  partial  differential  model  equations. 

Analysis  of  the  TSA  system  with  electrothermal  desorption  and  in-vessel  condensation  shows  that  the 
complete  cycle  can  be  broken  down  into  three  consecutive  steps: 

adsorption, 

desorption  without  condensation, 
desorption  accompanied  with  condensation. 

Although  these  three  stages  can  be  described  by  the  same  set  of  model  equations,  the  boundary 
conditions  corresponding  to  the  three  stages  are  essentially  different,  so  three  different  Femlab  models 
need  to  be  built.  In  order  to  simulate  a  complete  TSA  cycle,  these  3  models  have  to  be  run  in 
consequence  with  automatically  switching  from  one  model  to  the  next,  when  predefined  conditions  are 
reached.  The  most  convenient  way  to  achieve  that  is  to  build  a  combined  Femlab/Matlab  model  for  the 
complete  cycle.  Naturally,  separate  models  have  to  be  built  for  the  1 -cartridge  and  2-cartridges 
adsorbers. 

Taking  all  this  into  consideration,  these  are  the  main  phases  of  the  second  year  project: 

1.  Building  Femlab  models  for  adsorption,  desorption  and  desorption  accompanied  with 
condensation  for  1 -cartridge  adsorber; 

2.  Building  Femlab  models  for  adsorption,  desorption  and  desorption  accompanied  with 
condensation  for  2-cartridges  adsorber; 

3.  Building  Femlab/Matlab  models  for  integral  TSA  cycles  for  1 -cartridge  and  2-cartridges 
adsorbers; 

4.  Simulation  of  adsorption,  desorption  and  desorption-condensation  in  1 -cartridge  and  2- 
cartridges  adsorbers; 

5.  Simulation  of  complete  TSA  cycles  for  1 -cartridge  and  2-cartridges  adsorbers; 

6.  Investigation  of  the  influence  of  the  main  operational  parameters  of  the  TSA  process  and 
optimization  bases  on  it. 


2.  The  FEMLAB  models 

Six  Femlab  models  were  built: 

•  Model  A  l  -  model  for  adsorption  in  the  1-cartrdge  adsorber 

•  Model  D  l  -  model  for  desorption  without  condensation  in  the  1 -cartridge  adsorber 

•  Model  DC  l  -  model  for  desorption  with  condensation  in  the  1 -cartridge  adsorber 

•  Model_A_2  -  model  for  adsorption  in  the  2-cartrdges  adsorber 

•  Model_D_2  -  model  for  desorption  without  condensation  in  the  2-cartridges  adsorber 

•  Model_DC_2  -  model  for  desorption  with  condensation  in  the  2-cartridges  adsorber. 
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Model  assumptions: 

The  following  assumptions  were  used  in  setting  up  the  models: 

•  The  adsorbent  beds  are  treated  as  homogeneous,  with  uniform  adsorbent  porosity  and  density. 

•  The  mass  and  heat  transfer  resistances  on  the  particle  scale  are  neglected. 

•  The  fluid  phase  is  treated  as  an  ideal  gas  mixture  of  the  inert  and  the  adsorbate. 

•  All  physical  parameters  and  coefficients  are  considered  as  constants. 

•  The  electric  resistivity  of  the  ACFC  adsorbent  is  temperature  dependent.  Linear  temperature 
dependence  is  assumed,  based  on  experimental  results  reported  in  Ref  [5]. 

•  The  electric  power  during  electrothermal  desorption  is  supplied  under  constant  voltage 
conditions. 

•  The  condensation  at  the  adsorber  wall  is  dropwise.  This  assumption  is  based  on  experimental 
observations  reported  in  Ref  [5]. 

•  The  volume  of  the  condensed  liquid  is  neglected,  i.e.,  it  is  assumed  that  the  liquid  drops  don’t 
influence  the  gas  flow. 

•  The  heat  resistance  and  heat  capacity  of  the  adsorber  wall  are  neglected,  so  that  the  wall 
temperature  is  equal  to  the  temperature  of  the  environment. 

•  Initially,  the  adsorbate  concentrations  in  both  phases  are  in  equilibrium  and  uniform 
throughout  the  adsorbent  bed.  The  temperatures  of  both  phases  are  initially  equal  and  uniform 
throughout  the  adsorbent  bed. 


The  main  differences  from  the  models  given  in  Ref  [1 0] 

Owing  to  the  potentials  of  the  Femlab  software  to  model  rather  complex  systems,  some  of  the 
assumptions  used  during  the  first  phase  of  the  project  (Grant  No.  FA8655-03-1-3010  [10])  have  been 
omitted  it  the  models  presented  in  this  report.  These  are  the  basic  differences  between  the  current 
Femlab  models  and  the  models  developed  during  the  first  year  of  the  project: 

-  Two-dimensional  models  are  used  instead  of  one-dimensional  (the  gradients  both  in  the  axial 
and  in  the  radial  direction  are  taken  into  account). 

-  The  momentum  transport  equations  and  the  continuity  equations  for  the  inert  gas  are 
incorporated  in  the  models  and  solved  simultaneously  with  the  mass  and  energy  balance  equations. 

-  Instead  of  assuming  perfect  mixing  in  the  inlet  and  outlet  tubes,  the  velocity  field  is 
calculated  using  the  model  equations. 

-  The  gas  pressure  and  density  in  the  system  are  not  assumed  constant,  nor  uniform.  They  are 
calculated  from  the  model  equations. 

-  Condensation  at  the  adsorber  wall  is  taken  into  account. 


2.1.  Model  building  in  Femlab 

2.1.1.  Definition  of  the  system  geometry 

When  building  a  model  in  Femlab,  it  is  necessary  to  first  define  the  space  dimension.  Axially 
symmetrical  2D  space  was  used  for  both  adsorber  types. 

The  next  step  is  to  define  the  system  geometry.  This  is  done  by  direct  plotting  in  the  Femlab  window. 
The  system  geometry  is  shown  in  Figure  5:  Figure  5a  for  the  1 -cartridge  adsorber  and  Figure  5b  for 
the  2-cartridges  adsorber.  Figure  5  represents  the  halfs  of  the  axial  cross-sections  of  the  adsorbers 
(axial  symmetry  2D  models  are  used).  The  main  dimensions  of  the  adsorbers  are  also  shown. 
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(a) 


(b) 


Figure  5.  Definition  of  the  model  geometry,  with  the 
main  dimensions:  (a)  for  the  1-eartridge  and  (b)  for  the 
2-eartridge  adsorber 


2.1.2.  Definition  of  the  model  equations 

The  next  step  is  to  define  the  model  equations,  by 
choosing  the  appropriate  application  modes.  Figure  6 
shows  the  Model  Navigator  window  which  is  used  for 
defining  the  space  dimension  and  selecting  the 
application  modes. 

The  following  application  modes  were  used  in  our 
Femlab  models: 

1 .  Non-iso  thermal  flow  -  for  defining  the  gas  flow 
in  the  central  and  annular  tubes; 

2.  Brinkman  Equations  -  for  defining  the 
momentum  balance  for  the  packed  adsorbent 
bed(s); 

3.  Convection  and  Diffusion  -  for  defining  the 
mass  balance  for  the  gas  phase; 

4.  Diffusion  -  for  defining  the  mass  balance  for 
the  solid  phase; 

5.  Convection  and  Conduction  -  for  defining  the 
heat  balance  for  the  gas  phase; 


6.  Heat  transfer  by  Conduction  -  for  defining  the  heat  balance  for  the  solid  phase; 

7.  Conductive  Media  DC  -  for  defining  heat  generation  by  Joule  effect. 


Although  the  Coductive  Media  DC 
application  mode  is  not  needed  in  the 
adsorption  models  (there  is  no  electrical 
heating  during  adsorption),  all  7  application 
modes  were  used  in  all  models,  in  order  to 
enable  their  integration  into  the  models  of 
complete  TSA  cycles. 

The  geometry  of  the  adsorbers  assumes 
definition  of  their  parts,  i.e.  subdomains. 

Three  subdomains  are  defined  for  the  1- 
cartridge  adsorber  (see  Fig.  5a): 

Rl-  the  central  tube 
R2-  the  adsorbent  bed 
COl-  the  annular  space  around 
the  adsorbent  bed. 

On  the  other  hand,  5  subdomains  are  defined 
for  the  2-cartridges  adsorber  (Fig.  5b): 

Rl  -  the  central  tube  of  the 
bottom  cartridge 
R2  -  the  adsorbent  bed  of  the  bottom  cartridge 
R3  -  the  central  tube  of  the  top  cartridge 
R4  -  the  adsorbent  bed  of  the  top  cartridge 

COl  -  the  annular  space  between  the  cartridges  and  the  outer  adsorber  wall 


Figure  6.  The  Model  Navigator  window 


By  selecting  the  application  modes  one  by  one  in  the  Multiphysics  menu,  it  is  possible  to  access  the 
underlying  PDFs  and  define  their  coefficients  and  initial  conditions,  using  the  Subdomain  Settings 


5 


Menu  (an  example  is  shown  in  Figure  7),  and  the  boundary  conditions,  using  the  Boundary  Settings 
Menu  (an  example  is  shown  in  Figure  8),  for  each  subdomain. 

Checking  “Active  in  this  domain”  activates  the  chosen  application  mode  in  the  corresponding 
subdomain.  An  overview  of  the  application  modes  and  the  subdomains  in  which  they  are  active  is 
given  in  Table  1. 


Table  1.  An  overview  of  the  application  modes  and  the  subdomains  in  which  they  are  active 


Application  mode 

1-catridge  adsorber 

2-catridge  adsorber 

R1 

R2 

COl 

R1 

R2 

R3 

R4 

COl 

Non-isothermal  flow 

+ 

+ 

+ 

+ 

+ 

Brinkman  Equations 

+ 

+ 

+ 

Convection  and  Diffusion 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Diffusion 

+ 

+ 

+ 

Convection  and  Conduction 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Heat  transfer  by  Conduction 

+ 

+ 

+ 

Conductive  Media  DC 

+ 

+ 

+ 

Owing  to  the  complexity  of  our  model,  with  a  number  of  coupled  nonlinear  differential  equations  and 
interlinked  variables,  the  default  Femlab  PDFs  had  to  be  modified.  The  Femlab  window  Subdomain 
Settings  -  Equation  System,  such  as  the  one  shown  in  Figure  9,  is  used  for  this  purpose. 


Figure  8.  Boundary  Settings  Menu 


Figure  7.  Subdomain  Settings  Menu 
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Figure  9.  Subdomain  Settings  -  Equation  System  Menu 
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Figure  10.  The  initial  mesh  for  the  finite-differenee  method 
for:  (a)  1-eartridge  and  (b)  2-eartridges  models 


30 

After  modifying  the  PDEs  and  defining  the 
constants  and  the  parameters  in  the  system,  the 
mesh  for  the  finite  element  method  is 
generated.  Figure  10  shows  the  initial  meshes 
for  our  models  (Fig  10a  for  the  1 -cartridge  and 
Fig  10b  for  the  2-cartridges  adsorber).  ■  “ 

15 

The  models  were  solved  using  the  Direct 
(UMFPACK)  solver. 

to 
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(b) 


2.2.  The  underlying  equations 

The  following  equations  underlay  the  Femlab  models  presented  in  the  previous  section: 

1)  Momentum  balances  and  continuity  equations  for  the  inert  gas: 

For  the  central  and  outer  annular  tubes: 


du 
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(  ( du  ^  Sp  dp 

-  Pa  r  — + —  +u  =0 

^  \  dr  dz )  ]  dr  dz 


For  the  adsorbent  bed(s): 


du  dv 


+  ^  -2r|Li— +  -r|Li  — +  —  =-  r  -u  +  —  +2|li- 


dz  dr 


1^.. , 


k  dr 


dv  du 


dr  dz 


+  -M  ^  +  ^  -2r|i—  =-r  -v  +  — 


1^,. , 


k  dz 


du  dv 


^  I  Sr  dz 


-\  Pa\  n - ^ - - - =0 


dr  dz 


The  following  notations  are  used  in  these  equations  :  u  and  v  -  gas  veloeities  in  the  radial  and  in  the 
axial  direetion,  respeetively,  p  -  pressure,  Pg  -  gas  phase  density,  t  -  time,  r  -  radial  eoordinate,  z  - 
axial  eoordinate,  ^-adsorbent  bed  permeability  and  p  -  dynamie  viseosity  of  the  earrier  gas  (inert). 

The  adsorbent  bed  permeability  was  ealeulated  using  the  following  equation: 

k  =  (7) 

Vp 

The  gas  density  was  ealeulated  as: 

where  Tg  and  C  are  the  temperature  and  adsorbate  eoneentration  in  the  gas  phase  and  Rg  the  universal 
gas  eonstant. 


2)  Heat  balances  for  the  gas  phase: 

For  the  central  tube(s): 


-  '■pg  (^g  +  SvC)m  -  rpg  +  Cp.C)v-^ 


For  the  adsorbent  bed(s): 


r( *,o(r,  -  r, )  -  p.c,, J.f » ^  +  -  rp, (c„  +  c„,C)u ^  -  rp. (£■„  +  c,,C)v 


dr  dz 


For  the  outer  annular  tube: 


dC  dC 


-rpgCpJ\u  —  +  v—yrp^{Cp^+c^,C)u^-rp^{Cp^  +Cp,C)v^ 
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3)  Adsorbate  balance  for  the  gas  phase: 

For  the  central  tube(s): 


V71 

r - h  V| 

dt 


it 

„  it 

f 

Pg  Sr 

P. 

0zj 

V 

8C  dC^ 
ru - h  rv 


dr 


dz 


For  the  adsorbent  bed(s): 


dC  ^ 

r - h  V' 

dt 


f 


dC  dC 


V 


For  the  outer  annular  tube: 

A 


dz 


=  r- 


k^a 


(C  -C)- 


ru- 


dC  dC 


dr 


-  +  rv- 


dz 


dc  ^ 

r - h  Vi 

dt 


f 


-r- 


0C 

f 

r 

dr  Pg  dz  ^ 

V 

dC  dC^ 

ru - h  rv - 

dr  dz ) 


(12) 


(13) 


(14) 


4)  Electric  current  balance  for  resistive  heating: 


y  p  dr  p  dz 


=  0 


(15) 


In  this  equation  U  is  the  eleetrie  potential  and  p  the  eleetrie  resistivity  of  the  adsorbent,  whieh  is 
temperature  dependent.  Linear  temperature  dependenee,  whieh  was  obtained  experimentally  for 
ACFC  adsorbent  [5],  was  used  in  our  ealeulations: 


p=Po(i+Kr-7’«)) 


(16) 


5)  Heat  balance  for  the  solid  phase  within  the  adsorbent  bed(s): 


v|  rDt^  ^  +  rD, 


dt' 


I 


dr 


dz 


=  ^  -  rh,a{T^  -  T^) 

dV  dt 


(17) 


where  F  is  the  solid  temperature  and  /  JL)  the  eleetrie  power  supply  per  unite  volume  of  the 
adsorbent  bed: 


dV 


1 


^  dur  fdU^ 


2^ 


y  dr  J 


+ 


ydz  j 


1 


Po(1  +  M7;-7’«)) 


^  dUV  (dU^ 


2^ 


ydr  j 


+ 


ydz  j 


(18) 


6)  Adsorbate  balance  for  the  solid  phase  within  the  adsorbent  bed(s): 


p,^  =  k^a{C-C*)  (19) 

dt 

where  q  is  the  adsorbate  eoneentration  in  the  solid  phase  and  C*  the  adsorbate  eoneentration  in  the  gas 
phase  in  equilibrium  with  the  solid  phase.  The  Dubinin-Radushkevieh  equation  was  used  to  deseribe 
the  adsorption  equilibrium: 


f 

IL  =  ILoexp 

V 


Pa)) 


(20) 


In  this  equation  W  is  the  volume  of  the  adsorbed  phase  per  unite  mass  of  adsorbent,  W^-  the  volume  of 
the  mieropores  per  unit  mass  of  adsorbent,  E-  the  adsorbate  energy  of  adsorption,  p^-\hQ  adsorbate 
saturation  pressure,  -  the  adsorbate  partial  pressure.  The  saturation  pressure  was  ealeulated  using 
the  Wagner  equation: 
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(21) 


^  VP^x  +  VPgX^  -^  +  VPcX^  +  VPoX^  " 

UJ 

1  1-^  J 

where  x  =  \-{TjT^) ,  and  are  the  eritieal  pressure  and  temperature  and  VP  a,  VPb,  VPc  and  VPd 
are  Wagner  eonstants. 

Based  on  equation  (20),  the  following  relation  between  the  equilibrium  adsorbate  eoneentration  in  the 
gas  phase  C*  and  its  eoneentration  in  the  solid  phase  q  and  temperature  was  obtained: 


C 


pQXpl 


-In 


(22) 


(Mj  and  pj  are  the  adsorbate  molar  mass  and  density,  respeetively). 


Boundary  conditions: 

Equations  (1-19)  are  eommon  for  all  three  proeesses  (adsorption,  eleetrothermal  desorption  and 
eleetrothermal  desorption  with  eondensation)  and  for  both  adsorber  types  (1-eartridge  and  2-eartridges 
adsorbers),  and  eonsequently  for  all  six  Femlab  models.  Nevertheless,  some  of  the  boundary 
eonditions  for  these  6  models  are  different. 


Bl)  Common  boundary  conditions  for  all  models  for  the  1-cartridge  adsorber 

G 


z  =  0,  rG(0,q):  w=0,  v  = 


T  =T  C  =  C 

2  ^  g  gin  ’  ^  ’ 

'I  Tipg 


z  =  H  +  hi,  rG(0,rj):  u  =  0,  v  =  0,  -Dj’f]  = 


™  dz 


dz 


(23) 

(24) 


’  =  ri,  z&ih^,\+H):  P^^pI 

r  .  dT  ^  ^ 


Jit 


y  IzY  Qr 

f  .  QC  ^  ^ 

-D\.  — +  mC 


V 


\h  Qr 


+  h,fl-E,)(f-T^) 


Jb 


V 


\it 


dr 


Jit 


-D\—  +  uC 

’nr\b 


+  ^„,(1-8,)(C  -C), 


Jb 


-D':^  =  h,,{T^-TA,  J  =  0 


(25) 


dT„ 


dC 


z  =  h„  rs(p,rA:  u  =  0,  v  =  0,  -  D^fl  =  0,  -  D\^  =  0,U  =  U„  -  D^’;  =  0 


dP 


dz 


mz\b 


dz 


z  =  hi+H,  rs(p,r2):  m  =  0,  v  =  0, 


sn  , 


dz 


—  =  0,U  =  0,D^;^  =  0 
dz  dz 


dP  r)C^ 

zg(0, /![),  r  =  Tj  ;  u  =  0,  v  =  0,  -D^f\  —^  =  0,  =  0 


'it  dz 


mr\it 


dz 


(26) 

(27) 

(28) 


r  =  r2,  z&{h,h  +  H)\  p^=p^,,  v|^=v|^^, 


b  Qp 


ot 

\  r 


Jb 


-Atl^+P,(c,,+c  r)wr 


V 


\ot  Qp 


+  h^^(\-E,)(f-T\  (29) 


-D^]—  +  uC 

tttlb  Qj. 


\  f 
Jb  V 


■D^r\  ,  — +  mC 


+  42(1-s,)(C  -C),  J  =  0,  -D:;^  =  h^2iT,-TA 


dT 


lot 


dr 
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z  =  h  +  H,  re(r2,r4):  pl^=Pa 
Specific  boundary  conditions: 


v  =  0, 

-Di'fl 

=  0, 

-Z)„,  —  =  0 

\ot 

dz 

mz\ot 

\ot  ~  Pa  ’ 

dT^ 

dz 

=  0, 

(30) 

(31) 


For  the  1 -cartridge  adsorber,  the  basic  difference  between  the  three  models  is  in  the  boundary 
conditions  at  the  adsorber  wall  (for  r=rf)\ 

-  For  Model  A  l  and  Model  D  l  (adsorption  and  desorption  without  condensation)  they 
are  defined  as:  zero  velocity  at  the  wall,  zero  adsorbate  flux  through  the  wall  and  the  heat  flux  through 
the  wall  defined  by  the  heat  losses: 


r  =  r„  zs{h,h  +  H):  u  =  0,  v  =  0,  -DJ^—  =  0,  -  -TJ 


(32a) 


become: 


-  For  Model  DC  l  (desorption  with  condensation)  the  boundary  conditions  at  the  wall 


r  =  r 


4? 


zs{h,h  +  H):  u  =  0,v  =  0,C  =  C,JTJ,-Dt^ 


dT^ 

ot  Qr 


,  dC 

-^i-^cond)  +  KiT,-r)i'i2h) 


meaning  that,  when  the  conditions  for  condensation  at  the  wall  are  fulfilled,  the  concentration  of  the 
gas  phase  at  the  wall  becomes  equal  to  the  saturation  concentration  of  the  vapor  corresponding  to  the 
wall  temperature  7^,  Csat{Ty^>  At  the  same  time,  in  the  boundary  condition  for  the  heat  balance,  the 
heat  of  condensation  term  is  be  included  {{-IsHcond)  is  the  molar  heat  of  condensation).  In  our  model, 
the  heat  resistance  and  heat  capacity  of  the  wall  are  neglected  and  the  wall  temperature  is  assumed  to 
be  equal  to  the  ambient  temperature  Ta. 

B2)  Common  boundary  conditions  for  all  models  for  the  2-cartridges  adsorber 


z  =  0,  rG(0,q):  i/  =  0,  ^  =  — - 


T  =T  C=C 


z  =  h^+H,  re(0,ri):  u  =  Q,  v  =  0, 


dT„ 


dC 


z  =  h,+H  +  h„  re(0,r,):  u  =  0,  v  =  0,  -Z)f|  ^  =  0,  -DJ—  =  0 


'*  dz 


dT„ 


mzUt 

dC 


ze(0,/2i),  r  =  r^\  m  =  0,  v  =  0,  -Z)*^  ^  =  0,  -Z)„^|  ^  =  0 

dz  dz 

z^{h„H),  r  =  r,:  p|,=Z^|,,  m|.=“L’ 


r 


dT 

tr  I ..  Pje 


gy^pg  g 

\  J 

(  ,  sc  \  ( 

-Dj,  —  +  =  - 


it  Qp 


\  f 

it  V 


S7 


Jb 


(33) 

(34) 

(35) 

(36) 


n  \  ^  r 


+  k,,(l-s,)(C  -C), 


(37) 


=  J  =  0 

dr 


Jb 


dT 

z  =  h„  re(r„r2):  «  =  0,  v  =  0, -Z)^^[ -^  =  0, -Z)„,|^  ^  =  0,  Z/ =  t/o, -Z)“  ^  =  0 


dC 


dT, 

dz 


(38) 
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=  h,+H,  r^{r„r,)-.  «  =  0,  v  =  0, ^  =  0, ^  =  0,  C/ =  0, ^  =  0  (39) 


g(ZZ  +  /zi +/?2,2ZZ  +  /zi +/?2),  r  =  r^:  p\.,=p\^,,  “I,,  =«L’  H,-,  =Hi’ 


1,2  f'g'-  pg  '  g 


— +  mCJ  -|^-i)„„.|^  — +  mCJ^+A:„i(1-8;,)(C  -C), 
-D’:^^K,{T^-T,),  J  =  0 


=  h,+h,+H,  rs{r„r,):  «  =  0,  v  =  0, -Z)^^  ^  =  0, -Z)„, ,  ^  =  0,  Z/ =  Z/^, -Z)“  ^  =  0(41) 


9C_^rr_rr  r.hs  dT, 


M  ^_n  _n  I  ^-n  ..-n  _  ^ 


z  =  2H  +  h,+h„  rs(r„r,):  „  =  0,  v  =  0, -Z),^^|  ^  =  0, -Z)„,  —  =  0,  Z/ =  0, -Z),f  ^  =  0  (42) 


zg(0,/!i),  r  =  r2:  m  =  0,  v  =  0,  -Dj’fl  ^  =  0,  -D\  —  =  0 


z^{h„k+H),  r^r,:  p^^p^^, 


-D’:^\^+p(c^^+cC)uT  =  ^+p  (c  +c  nwr  +/^^^(i-e^)(r^-r), 


+K2i^-^b)ic  -c),  J^o, 

dr  Jr  \  or  j. 


K^-KiiT.-T,) 


zs(h,+H,h,+h,+H),  r  =  r,:  u  =  0,  v  =  0,  -D^^l  ^  =  0,  -D^^—  =  0 


zg(/?i  +  /?2+ZZ,/?i+/22  +  2ZZ),  r  =  r^-.  Z?|^=p[,,  mL=mL,, 

^  dT  ^  r  dT  ^ 

=  -^^l^+p.(s.  +  SvO«^J  +/z,2(i-£,)(7;-rj, 


lot  Qr  ^  P'’  '  g 


-D„^—  +  uC  =  -D„^—  +  uC  +k„,{\-z,){0-C),  J  =  0, 


’nr\b 


-d':^=k,{t^-t^) 


z  =  2H  +  hi  +  h2,  rs{r2,r^):  u  =  0,  v  =  0,  -  Dj'f  — ^  =  0,  - =  0 


zg(0,/!i),  re{r2,r^):  u  =  0,v  =  0,  - =  0,- =  0  (48) 


Specific  boundary  conditions: 

For  the  2-caartridges  adsorber,  the  differences  between  the  models  corresponding  to  adsorption  and 
desorption  are  a  result  of  different  gas  flow  directions  in  these  two  steps.  Accordingly,  the  boundary 
conditions  for  the  adsorber  inlets  and  outlets  are  different.  On  the  other  hand  the  boundary  conditions 
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at  the  adsorber  wall  differ  for  the  models  without  and  with  condensation  (in  an  analogous  way  as  for 
the  1 -cartridge  models).  The  corresponding  boundary  conditions  are: 

-  For  Model  A  2: 


z  =  0,  re(r2,r3):  ^  =0, 

^ot  dz 

-A.Lf  =  o, 

II 

O 

II 

o 

(49a) 

z-\+h2+2H,  rE(0,q):  p-Pa-> 

-D  1  — =0 

(50a) 

-  For  Model_D_2  and  Model_DC_2: 

dr 


z  =  0,  = 

ot  dz 


^  =  0 

mz\nt  P  Pa 


(49b) 


z  =  h^+h2+2H,  re(0,r^):  u  =  0, 


G 

2  ’ 

^Pg 


T  =T. 

g  gin 


C  =  G 


(50b) 


-For  Model  A  2  and  Model  D  2: 


zs(h„2H  +  h,+h,),  r  =  r,:  u  =  0,v  =  0,- =  =  -TJ 


(51a) 


-For Model  DC  2: 


z  &  (hi,2H  +  hi  +  h2),  r  =  r^\  u  =0,v  =  0,C  =  C^^,{TJ, 


8T„  .  dC 

tr  ~^'nr\ot  cond  )  +  Kg  (^g  “  ) 


(51b) 


In  equations  (1-51),  G  is  the  molar  flow-rate  of  the  gas  inert,  Uq  -  the  supply  voltage  and  Sb  -  bed 
porosity.  The  subseript  in  denotes  the  inlet,  o  the  outlet,  it  the  inner  (eentral)  tube,  ot  the  outer 
(annular)  tube  and  a  the  ambient.  The  meanings  of  the  geometrie  values  {r\  to  r4,  H,  h\  and  are 
defined  in  Figure  5.  A  eomplete  list  of  the  notations,  ineluding  the  physieal  parameters  and 
eoeffieients  used  in  the  model  equations,  ean  be  found  in  the  nomenelature,  at  the  end  of  this  report. 

Although  the  heat  and  mass  transfer  eoeffieients,  dispersion  eoeffieients,  ete.,  ean  be  defined  as 
funetions  of  the  dependent  variables  (veloeities,  eoneentrations,  temperatures),  in  our  eurrent  models 
eonstant  values  were  used  in  order  to  reduee  the  eonvergenee  problems.  The  adsorption  isotherm 
relation  {q=i{C,  T^)  and  the  temperature  dependenee  of  the  eleetrie  resistivity  were  ineluded  in  the 
models  (equations  (22)  and  (16),  respeetively). 
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3.  Models  for  complete  TSA  cycles 

In  order  to  model  a  complete  TSA  cycle,  it  is  necessary  to  switch  from  one  Femlab  model  to  another, 
when  certain  conditions  are  fulfilled. 

In  order  to  reduce  the  time  needed  to  achieve  the  periodic  quasi-steady  state,  the  TSA  cycle  was 
started  with  the  desorption  step.  The  initial  state  of  the  system  corresponds  to  equilibrium  at  ambient 
temperature.  The  initial  concentration  in  the  gas  phase  is  equal  to  the  feed  concentration  during 
adsorption: 

?<0,Vr,Vz:  T^=T^=T^.^=T^^=T^=T^,C  =  C,„=C*  =^{q)  =  C^,u  =  v  =  Q,  p  = 


Three  switching  moments  are  defined: 


Time=21  [s]  Surface:Concentration  Subdomain  marker: concentration 


20 


0  ■  • 
-10 


condensation 

started 


=  0.0104462 


r[cm] 


1.  Switching  from  ModelD  to  ModelDC. 

This  switch  is  performed  at  the  moment  when 
the  gas  concentration  at  the  adsorber  wall  (for 
r=rf)  becomes  equal  to  the  saturation 
concentration  corresponding  to  the  wall 
temperature  Csat,  resulting  with  start  of 
condensation.  The  end  of  condensation  can  be 
detected  from  the  change  of  the  sign  of  the 
adsorbate  flux  at  the  adsorber  wall 

{- becomes  <0).  In  Figure 

11  we  show  the  distribution  of  the  gas 
concentrations  in  subdomain  COl  at  the  first 
switching  moment,  for  the  1-cartidge  adsorber. 

It  can  be  seen  that  the  concentration  is  nearly 
uniform  at  the  wall  surface,  meaning  that  the 
concentration  practically  starts  on  the  whole 

surface  simultaneously.  A  similar  result  is  obtained  for  the  2-cartridges  adsorber. 

2.  Switching  from  Model  DC  to  Model_A. 

This  switch  defines  the  end  of  the  desorption  and 
the  start  of  the  adsorption  part  of  the  cycle.  It  is 
performed  when  the  maximal  temperature  in  the 
adsorbent  bed  reaches  certain  predefined  value 
Tsw,  which  should  not  be  exceeded.  The 
temperature  distribution  in  the  solid  phase  at  the 
moment  when  the  maximal  temperature  reaches 
473.15  K  is  shown  in  Figure  12,  for  the  1-cartridge 
adsorber.  It  can  be  seen  that  significant 
temperature  gradients  are  present  both  in  the  axial 
and  in  the  radial  directions  in  the  adsorbent  bed. 


iiri  =  0.0104251 


10.0104 


0.0104 


■  ^0.0104 

10.0104 
0.0104 


Figure  11.  Maximal  and  minimal  concentration  at  the 
adsorber  wall  in  the  C  switching  moment  (1-catridge) 
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3.  Switching  from  Model_A  to  Model  D. 

This  switch  corresponds  to  the  end  of  adsorption 
and  start  of  the  new  desorption  step.  It  is 
performed  automatically  when  the  outlet 
concentration  reaches  a  certain  predefined  breakthrough  value  Cbreab 
concentration  is  obtained  as  the  average  value  across  the  outlet: 

r.  —  To 

3  2  Vr, 


r[cm] 

Figure  12.  Maximal  and  minimal  temperature  in  the 
adsorbent  bed  in  the  2^^  switching  moment  (1-catridge) 

In  principle,  the  outlet 
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Time=2884  [s]  Surface:Concentration  Subdomain  marker: concentrationx  iq'** 

From  Figure  13,  which  shows  the  concentrations 
in  the  annular  tube  of  the  1-catridge  adsorber,  at 
the  moment  of  breakthrough,  it  can  be  seen  that 
the  concentration  is  practically  uniform  across  the 
outlet.  The  same  result  is  obtained  for  the  2- 
cartridges  adsorber. 

Although  in  principle  different  Femlab  models  can 
be  run  consecutively,  automatic  checking  whether 
certain  conditions  are  met  and  switching  from  one 
model  to  another  based  on  it,  is  not  possible  in 
Femlab.  For  that  reason,  building  of  the  integral 
models  for  complete  TSA  cycles  was  done  by 
using  combination  of  Femlab  and  Matlab. 

The  procedure  for  building  these  models  can 
shortly  be  given  by  the  following  steps: 

1.  Building  and  running  the  Femlab  models  for  separate  stages  of  the  TSA  cycle  (Model  A  l, 
Model  D  l  and  Model  DC  l  for  1 -cartridge  and  Model_A_2,  Model_D_2  and  Model_DC_2 
for  2-cartridges  adsorber).  The  same  application  modes  have  to  be  used  in  all  models.  The 
same  sequence  of  the  application  modes  and  of  the  subdomains  has  to  be  used  in  all  models 
corresponding  to  the  same  type  of  adsorber. 

2.  Saving  the  Femlab  models  as  Matlab  programs  (m-files). 

3.  Integrating  all  three  m-files  into  one  Matlab  program  (some  program  segments,  like  geometry 
definition  and  mesh  initialization,  are  needed  only  once). 

4.  Adding  the  necessary  Matlab  commands: 

a.  For  checking  whether  certain  conditions  are  met; 

b.  For  defining  the  new  initial  state  as  the  solution  of  the  previous  step; 

c.  For  switching  from  one  to  another  model; 

d.  For  repeating  the  complete  cycle  the  needed  number  of  times. 

5.  Plotting  the  results  in  Matlab  or  exporting  them  as  a  fern  structure  that  can  be  imported  into 
Femlab  for  postprocessing. 

The  structure  of  the  integral  TSA  models  is  given  in  Figure  14.  This  structure  is  common  for  both 
configurations  of  the  adsorber  (1 -cartridge  and  2-cartridges),  so  common  names  Model  D,  Model  DC 
and  Model  A  are  used. 
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Figure  13.  Maximal  and  minimal  concentration  at  the  outlet 
in  the  3^^  switching  moment  (1-catridge) 
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Figure  14.  Structure  of  the  complete  TSA  model 
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4.  Simulation  results 


The  Femlab  and  Matlab/Femlab  models  were  used  for  simulation  of  adsorption,  eleetrothermal 
desorption  without  and  with  eondensation  and  of  eomplete  TSA  eyeles. 

Simulation  of  adsorption  was  performed  for  initially  elean  adsorbent  {qp=0)  at  ambient  temperature, 
and  for  eonstant  feed  eoneentration  and  flow-rate. 

Simulation  of  desorption  was  performed  for  initially  uniformly  saturated  adsorbent  at  ambient 
temperature,  for  eonstant  voltage  supply  and  eonstant  flow-rate  inert  feed  stream  (C/„=0).  The 
simulation  results  presented  in  this  report  eorrespond  to  a  eomplete  desorption  (regeneration)  step, 
ineluding  desorption  before  the  start  of  eondensation  and  desorption  aeeompanied  with  eondensation. 

The  simulations  were  performed  for  a  laboratory  seale  adsorber  of  the  following  dimensions: 

-  ri=0.95  em 

-  r2=1.5  em 

-  r^=2  em 

-  r4=3.5  em 

-  if=30  em 

-  h\=2  em 

-  h2=2  em 

and  for  the  following  system: 

-  Adsorbate:  methyl  ethyl  ketone  (MEK) 

-  Adsorbent:  Ameriean  Kynol  ACC-5092-20  (aetivated  earbon  fiber  eloth  material) 

-  Carrier  gas:  nitrogen. 

The  bed  permeability  was  ealeulated  from  the  experimentally  measured  pressure  drop  of  the  ACFC 
bed  given  in  Ref  [5].  The  numerieal  values  of  most  of  the  other  model  parameters  are  also  based  on 
Ref  [5].  A  eomplete  list  of  the  model  parameters  used  in  these  simulations  is  given  in  Table  2. 

The  simulations  were  performed  with  the  initial  meshes  generated  in  Femlab,  and  shown  in  Figures 
10a  and  10b. 

Femlab  offers  a  large  number  of  options  for  graphieal  representation  of  the  simulated  results.  The 
following  results  were  ehosen  for  presentation  in  this  report: 

Veloeity  field  in  the  adsorber; 

Gas  and  solid  eoneentration  fields  in  the  adsorber; 

Gas  and  solid  temperature  fields  in  the  adsorber. 

These  results  are  shown  in  different  forms:  as  2-D  surfaee  plots  and  as  time,  radial  and  axial  profiles, 
in  1-D  diagrams. 

For  desorption,  the  eleetrie  resistivity  field  is  also  shown. 
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Table  2.  The  model  parameters  used  for  simulation 


Parameter 

Notation  and  value  | 

Adsorption 

Desorption  | 

Adsorber  dimensions 

ri=0.95  em  r2=1.5em 

r3=2  em  r4=3.5  em 

iT=30  em  h\^  /z2=2  em 

Bed  porosity 

8^=0.72 

Speeifie  surfaee  area 

a=13.65  em^/em^ 

Ambient  pressure 

Patm=101325  Pa 

Inlet  adsorbate  eoneentration  in  the  gas  phase 

Q„=0.001  mol/mol 

mol/mol 

Initial  adsorbate  eoneentration  in  the  gas  phase 

C^=0  mol/mol 

C^=0.001  mol/mol 

Initial  adsorbate  eoneentration  in  the  solid  phase 

qp^O  mol/g 

^^=0.00465  mol/g 

Saturation  eoneentration 

C^^^=0.01  mol/mol  | 

Eleetrie  voltage  supply  Uq 

f/o=0  V 

f/o=llV 

Inert  gas  flow-rate 

G=0.1  mol/s 

G=0.01  mol/s 

Switehing  temperature  from  desorption  to  adsorption 

r,^=463.15  K 

Breakthrough  eoneentration  for  switehing  from  adsorption  to  desorption 

Qreayt^O-0005  mol/mol 

Inlet  gas  temperature 

re,„=293.15  K 

Ambient  temperature 

r^=293.15  K 

Wall  temperature 

r^=293.15  K 

Mass  transfer  eoeffieient  in  the  bed 

^^=0.0003  mol/(em^s) 

Bed  to  eentral  tube  mass  transfer  eoeffieient 

^^7=0.0003  mol/(em^s) 

Bed  to  annular  tube  mass  transfer  eoeffieient 

^^2"'0.0003  mol/(em^s) 

Heat  transfer  eoeffieient  in  the  bed 

**=0.0055  W/Ccm^s) 

Bed  to  eentral  tube  heat  transfer  eoeffieient 

*,/=0.001  W/(cm^s) 

Bed  to  eentral  tube  heat  transfer  eoeffieient 

*j2=0.001  W/(cm^s) 

Heat  transfer  eoeffieient  defining  heat  losses 

*,„e=5xl0-^  W/(cm^K) 

Radial  dispersion  eoeffieient  of  the  gas  phase  in  the  eentral  and  annular  tubes 

-DmJ  it=£>n,rl  ot=l  mol/(cms) 

Axial  dispersion  eoeffieient  of  the  gas  phase  in  the  eentral  and  annular  tubes 

-^mzl  it  -^mzl  ot  1  mol/(ems) 

Radial  heat  diffusivity  of  the  gas  phase  in  the  eentral  and  annular  tubes 

Ar‘'®l  it=Ar''®l  ot=l  W/(Kcm) 

Axial  heat  diffusivity  of  the  gas  phase  in  the  eentral  and  annular  tubes 

Az"®!  ot=l  W/(Kcm) 

Radial  and  axial  heat  diffusivity  of  the  solid  phase 

A* -l-2xl0‘^  W/(cmK) 

Molar  heat  of  adsorption 

(-A/4j,)=66200  J/mol 

Molar  heat  of  eondensation 

(-AA,„„4=31230  J/mol 

Adsorbent  bed  density 

p/,=0.221  g/em^ 

Adsorbate  density 

p^=0.81  g/em^ 

Adsorbate  molar  mass 

M^=72.107  g/mol 

Inert  gas  molar  mass  (nitrogen) 

M5=28.02  g/mol 

Speeifie  heat  eapaeity  of  liquid  adsorbate 

c„/=157.9  J/(molK) 

Speeifie  heat  eapaeity  of  the  inert  gas 

c^e=29.13  J/(molK) 

Heat  eapaeity  of  the  solid  phase 

c„,=0.71  J/(gK) 

Speeifie  heat  eapaeity  of  gaseous  adsorbate 

c^,=100.9  J/(molK) 

Dynamie  viseosity  of  the  inert  gas 

p=L78xlO'"^  g/em/s 

Inert  gas  thermal  eonduetivity 

*=2.6x1  O''*  W/cm/K 

Eleetrie  resistivity  at  referent  temperature 

po=0.202  flem 

Temperature  eoeffieient  of  the  bed  eleetrieal  resistivity 

*=-2x10'^  1/K 

Critieal  temperature  of  the  inert  gas  (nitrogen) 

T,b=\262  K 

Critieal  pressure  of  the  adsorbate 

/7<,zi=4.26x10+‘’  Pa 

Critieal  pressure  of  the  inert  gas  (nitrogen) 

^5,5=3.39x10^®  Pa 

Wagner  eonstants  for  the  adsorbate 

FP^=-7.71476,  VPb=\.1\06\ 
FPc^-3.6877,  VP .15X69 

Total  volume  of  mieropores  (D-R  equation) 

lJo=0.748  cm^/g 

Boltzmann’s  eonstant 

4=1.38048x10'“  J/K 

Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

£'=14.43x10'"^  J/mol 
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Two  additional  parameters,  very  important  for  the  eleetrothermal  desorption  proeess,  were  ealeulated: 
the  mass  of  the  eondensed  liquid  and  the  used  eleetrie  energy. 


Calculation  of  the  mass  of  the  condensed  liquid: 

The  flux  of  the  eondensed  liquid  is  ealeulated  from  the  gas  eoneentration  gradient  at  the  adsorber  wall: 


,  .  _  I  9C| 


(52) 


This  flux  is  a  fiinetion  of  the  axial  position  and  time.  The  eondensation  rate  is  ealeulated  by  integration 
of  the  flux  over  the  whole  surfaee  of  the  adsorber  wall,  and  is  a  fiinetion  of  time  only: 

Htot 

kond (0  =  27ir4  j  (z,  t)dz  (53) 

0 

The  total  amount  of  the  eondensed  liquid  is  obtained  by  integrating  the  eondensation  rate  over  time: 

'^cond 

Kond=  \Lcondd)dt  (54) 

0 

Htot  denotes  the  total  height  of  the  adsorber  and  Xcond  the  total  time  of  condensation. 


Calculation  of  the  used  electric  energy: 

It  is  assumed  that  the  used  eleetrie  energy  is  equal  to  the  Joule  heat  produeed  in  the  adsorbent  material. 
This  energy  per  unit  volume  and  unit  time  is  defined  by  equation  (18). 


The  used  eleetrie  power  is  obtained  by  integration  over  the  adsorbent  bed  volume: 


Qelit)- 
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f  f  r  ^ 

y  56/ ) 
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L  4 

1  dr  ) 

t  5z  J 
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liidrdz 


(55) 


and  the  total  used  energy  during  desorption,  by  integration  over  time: 


'^des 

Qei=jai(0dt  (56) 

0 

Tdes  is  the  total  desorption  time. 


4.1.  Simulation  results  for  the  1-cartridge  adsorber 

4.1. 1.  Adsorption  in  the  1-cartridge  adsorber 

The  simulation  of  adsorption  is  performed  for  7000  s.  After  that  time  the  eolumn  is  praetieally 
saturated  and  the  adsorption  is  finished.  The  results  presented  here  were  obtained  for  inlet 
eoneentration  0.001  mol/mol  and  inert  gas  flow  rate  0.1  mol/s. 

Velocity  distribution 

The  veloeity  field  in  the  1-eartridge  adsorber  is  presented  first.  In  Figure  15,  the  veloeity  distribution 
at  the  end  of  adsorption  (at  t=7000  s)  is  shown  in  the  form  of  a  2-D  surfaee  plot,  together  with  a  veetor 
representation  of  the  veloeity.  Constant  veloeity  field  was  obtained  nearly  instantaneously  (during  the 
first  3  seeonds).  The  veloeity  distribution  is  also  presented  in  1-D  diagrams:  the  time  profiles  at  the 
middle  eross-seetion  of  the  adsorber  (z=l6  em)  and  4  radial  positions  in  Figure  16,  the  radial  profiles 


19 


for  ^=7000  s  and  4  axial  positions  in  Figure  17  and  the  axial  profiles  for  ^=7000  s  and  4  radial  positions 
in  Figure  18.  Some  numerieal  noise  is  observed  in  the  simulated  gas  veloeities  in  the  eentral  tube. 
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Figure  15.  Velocity  distribution  of  the  gas 
phase  during  adsorption  in  the  1 -cartridge 
adsorber,  for  /=7000s 


Figure  16.  Time  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  at  z=  16  cm 


Figure  17.  Radial  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  for  /=7000  s 


Figure  18.  Axial  profiles  of  the  gas  velocity  during 
adsorption  in  the  1 -cartridge  adsorber,  for  t=7000  s 


Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  adsorption  in  the  1 -cartridge  adsorber  are  first 
presented  in  the  form  of  2-D  surface  plots,  in  Figures  19  and  20  respectively.  The  concentrations  at 
three  distinct  moments  (^=50,  200  and  7000  s)  are  shown,  in  order  to  illustrate  their  evolution  in  time. 
For  ^=7000  s  the  gas  concentration  in  the  adsorber  is  nearly  uniform  and  equal  to  the  feed 
concentration  (0.001  mol/mol),  which  means  that  the  adsorption  process  is  finished.  At  early  stages  of 
adsorption,  the  concentrations  in  the  adsorber  are  nonuniform,  both  in  the  axial  and  in  the  radial 
direction. 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  16  cm  in  Figure  21, 
the  radial  profiles  at  z=16  cm  in  Figure  22  and  the  axial  profiles  for  ^=200s  in  Figure  23.  Figures  21a, 
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22a  and  23a  correspond  to  the  gas  phase  and  Figures  21b,  22b  and  23b  to  the  solid  phase.  For  the 
model  parameters  used  in  our  simulations,  the  gas  concentration  radial  profiles  in  the  central  and  in  the 
annular  tube  are  practically  flat. 
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Figure  19.  Gas  concentration  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  /=50s,  (b)  /=200s  and 
(c)  /=7000s 
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Time=50  s  Surface:  Solid  concentration  Time=200  s  Surface:  Solid  concent  Time=7000  s  Surface:  Solid  concent 


Figure  20.  Solid  concentration  distribution  during  adsorption  in  the  1 -cartridge  adsorber:  (a)  /=50s,  (b)  /=200s  and 
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(c)  /=7000s 
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Figure  21.  Time  profiles  of  the  eoneentrations  during  adsorption  in  the  1-eartridge  adsorber,  at  z=16em 
(a)  gas  phase;  (b)  solid  phase 
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hgure  22.  Radial  profiles  of  the  eoneentrations  during  adsorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  23.  Axial  profiles  of  the  eoneentrations  during  adsorption  in  the  1-eartridge  adsorber,  for  /=200  s: 
(a)  gas  phase;  (b)  solid  phase 


22 


Temperature  distributions 


The  gas  and  solid  temperature  distributions  in  the  1-eartridge  adsorber  are  first  given  in  Figures  24  and 
25,  respeetively,  in  the  form  of  2-D  surfaee  plots.  The  simulated  gas  and  solid  temperatures  are  also 
shown  in  1-D  diagrams:  the  time  profiles  in  Figures  26a  and  26b,  the  radial  profiles  in  Figures  27a  and 
27b  and  the  axial  profiles  in  Figures  28a  and  28b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases.  As  the  temperature  ehanges  in  the  adsorber  are  due  only  to  the  heat 
generated  by  adsorption,  they  are  not  signifieant  (the  maximal  temperature  ehange  is  less  than  5  ®C. 
The  biggest  temperature  ehange  is  observed  during  the  first  200  s  of  the  adsorption  process. 
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Figure  24.  Gas  temperature  distribution  during  adsorption  in  the  1-eartridge  adsorber:  (a)  /=50s,  (b)  /=200s  and  (e)  /=7000s 
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Figure  25.  Solid  temperature  distribution  during  adsorption  in  the  1-eartridge  adsorber:  (a)  /=50s,  (b)  /=200s  and  (e)  /=7000s 
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Figure  26.  Time  profiles  of  the  temperatures  during  adsorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  26.  Radial  profiles  of  the  temperatures  during  adsorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 

(a)  (b) 


Figure  28.  Axial  profiles  of  the  temperatures  during  adsorption  in  the  1-eartridge  adsorber,  for  /=200  s 
(a)  gas  phase;  (b)  solid  phase 


24 


4.1.2.  Desorption  in  the  1-cartridge  adsorber 

The  complete  desorption  step,  before  and  after  the  start  of  condensation,  was  simulated.  This  was  done 
by  using  a  Matlab  model  obtained  by  combining  Model  D  l  and  Model  DC  l.  The  results  presented 
here  were  obtained  for  initial  concentration  in  the  solid  phase  0.00465  mol/g  (equilibrium 
concentration  for  C=0.001  mol/mol),  electric  voltage  11  V  and  inert  gas  flow-rate  0.01  mol/s.  The 
time  scale  used  for  this  simulation  was  0  to  200  s.  After  200  s  the  desorption  process  is  nearly 
finished. 

Velocity  distribution 

The  velocity  distribution  at  the  end  of  desorption  (at  ^=200  s)  is  shown  in  Figure  29,  in  the  form  of  a  2- 
D  surface  plot,  together  with  a  vector  representation  of  the  velocity.  The  velocity  distribution  is  also 
presented  in  1-D  diagrams:  the  time  profiles  in  Figure  30,  the  radial  profiles  in  Figure  31  and  the  axial 
profiles  in  Figure  32.  Some  (not  significant)  change  of  the  velocity  profile  with  time  is  observed,  as  a 
result  of  the  change  of  the  gas  temperature. 
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Figure  29.  Velocity  distribution 
during  desorption  in  the  1 -cartridge 
adsorber,  for  t=200  s 


t[s] 

Figure  30.  Time  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  at  z=  16  cm 


Figure  31.  Radial  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  for  ^200  s 


Figure  32.  Axial  profiles  of  the  gas  velocity  during 
desorption  in  the  1 -cartridge  adsorber,  for  /=  200  s 
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Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  desorption  in  the  1 -cartridge  adsorber,  presented 
in  the  form  of  2-D  surface  plots,  are  given  in  Figures  33  and  34  respectively.  The  concentrations  at 
three  distinct  moments  (^=24,  80  and  200  s)  are  shown.  ^=24  s  is  the  moment  when  condensation 
starts,  i.e.  when  the  program  switches  from  Model  D  l  to  Model  DC  l.  After  that  time  the 
concentration  at  the  adsorber  wall  (and  practically  in  the  whole  annular  tube)  becomes  constant  and 
equal  to  the  saturation  concentration  at  293.15  K  (0.01  mol/mol  for  MEK). 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  z=  16  cm  in  Figure  35, 
the  radial  profiles  at  z=16  cm  in  Figure  36  and  the  axial  profiles  for  ^  =  80  s  in  Figure  37.  Figures  35a, 
36a  and  37a  correspond  to  the  gas  phase  and  Figures  35b,  36b  and  37b  to  the  solid  phase. 
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Figure  33.  Gas  concentration  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  /=24s,  (b)  /=80s  and 
(c)  /=200s 
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Figure  34.  Solid  concentration  distribution  during  desorption  in  the  1 -cartridge  adsorber:  (a)  /=24s,  (b)  /=80s  and 
(c)  /=200s 
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Figure  35.  Time  profiles  of  the  eoneentrations  during  desorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  36.  Radial  profiles  of  the  eoneentration  during  desorption  in  the  1-eartridge  adsorber,  at  z=16em: 
(a)  gas  phase;  (b)  solid  phase 
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Figure  37.  Axial  profiles  of  the  eoneentrations  during  desorption  in  the  1-eartridge  adsorber,  for  /=80  s: 
(a)  gas  phase;  (b)  solid  phase 
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Condensation  flux  and  condensation  rate 

The  condensation  flux  is  calculated  using  equation  (52).  It  is  a  function  of  both  time  and  axial  position 
in  the  adsorber.  The  time  profiles  of  the  condensation  flux  are  shown  in  Figure  38,  and  the  axial 
profiles  in  Figure  39. 


Figure  38.  Time  profiles  of  the  eondensation  flux  in 
the  1-eartridge  adsorber 


These  figures  show  that  condensation  starts  at 
t=2A  s  and  the  condensation  flux  is  maximal  at 
approximately  ^=80s.  At  ^=200s,  when  the  power 
is  switched  of,  the  condensation  is  still  not 
completely  finished.  At  the  beginning  of  the 
process  the  condensation  flux  increases  from  the 
bottom  to  the  top  of  the  adsorber,  but  with  time 
the  axial  profile  changes  and  the  condensation  is 
most  intensive  near  the  bottom  of  the  adsorber. 

Figure  40  shows  the  time  profile  of  the 
condensation  rate,  calculated  by  integration  of 
the  condensation  flux  along  the  adsorber  wall 
surface  (equation  (55)). 


Figure  39.  Axial  profiles  of  the  eondensation  flux  in 
the  1-eartridge  adsorber 


TIME  (s) 


Figure  40.  Condensation  rate  in  the  1-eartridge  adsorber 


Temperature  distributions 

The  2_D  surface  plots  representing  the  gas  and  solid  temperature  distributions  in  the  1 -cartridge 
adsorber  are  given  in  Figures  41  and  42,  respectively.  The  simulated  gas  and  solid  temperatures  are 
also  shown  in  1-D  diagrams:  the  time  profiles  in  Figures  43a  and  43b,  the  radial  profiles  in  Figures 
44a  and  44b  and  the  axial  profiles  in  Figures  45a  and  45b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases  and  both  directions.  During  desorption,  the  temperatures  in  the 
adsorber  change  significantly,  owing  to  the  production  of  Joule  heat. 
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Figure  41.  Gas  temperature  distribution  during  desorption  in  the  1-eartridge  adsorber:  (a)  /=24s,  (b)  /=80s  and 
(e)  /=200s 
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Time=80  Surface:  Solid  temperature 
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Figure  42.  Solid  temperature  distribution  during  desorption  in  the  1-eartridge  adsorber:  (a)  /=24s,  (b)  /=80s  and  (e)  /=200s 
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Figure  43.  Time  profiles  of  the  temperatures  during  desorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  44.  Radial  profiles  of  the  temperatures  during  desorption  in  the  1-eartridge  adsorber,  atz=16em: 
(a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  45.  Axial  profiles  of  the  temperatures  during  desorption  in  the  1-eartridge  adsorber,  at  r=\25  em: 
(a)  gas  phase;  (b)  solid  phase 
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Distribution  of  the  electric  resistivity  of  the  ACFC  adsorbent  bed 

The  electric  resistivity  of  the  ACFC  is  temperature  dependent  (decreases  with  the  increase  of 
temperature).  As  the  temperature  in  the  adsorbent  bed  varies  considerably  both  in  time  and  space,  so 
does  the  electric  resistivity.  This  dependence  is  shown  in  Figure  46  (in  the  form  of  a  2-D  surface  plot) 
and  Figures  47,  48  and  49  (1-D  diagrams  showing  time,  radial  and  axial  distributions,  respectively). 

The  electric  power  transformed  into  Joule  heat  also  changes,  as  a  consequence  of  the  change  of  the 
electric  resistance.  For  constant  voltage  supply,  the  electric  power  increases  with  time,  as  shown  in 
Figure  50.  The  electric  power  is  calculated  using  equation  (55). 
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Figure  46.  The  eleetrie  resistivity  of  the  ACFC  adsorbent  distribution  in  the  1-eartridge  adsorber:  (a)  /=24s,  (b)  /=80s  and 
(e)  /=200s 


Figure  47.  Time  profiles  of  the  eleetrie  resistivity 
during  desorption  in  the  1-eartridge  adsorber,  at 
z=16em 


Figure  48.  Radial  profiles  of  the  eleetrie  resistivity 
during  desorption  in  the  1-eartridge  adsorber,  at 
z=16em 
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Figure  49.  Axial  profiles  of  the  eleetrie  resistivity  in 
the  1-eartridge  adsorber,  for  t=80  s 


Figure  50.  Eleetrie  power  supply  during  desorption 
in  the  1  -eartridge  adsorber 


4.1.3.  TSA  process  in  the  1-cartridge  adsorber 

The  complete  TSA  process  was  simulated  by 
using  the  integral  Matlab/Femlab  model  presented 
in  Figure  14.  These  simulations  take  long  times,  so 
they  were  limited  to  the  first  three  cycles.  The 
simulation  results  presented  here  were  obtained  for 
7^^=463.15  K  and  Cbreak=0.5  Cfn.  All  other 
simulation  parameters  are  the  same  as  in  the 
simulations  of  adsorption  and  desorption.  The  time 
and  space  distributions  of  the  velocity, 
concentrations  and  temperatures  in  the  system 
during  the  adsorption  and  desorption  half-cycles 
are  very  similar  to  the  results  presented  in  Sections 
4.1.1  and  4.1.2,  so  they  will  not  be  given  here. 
Instead,  we  are  presenting  the  change  of  the  outlet 
concentration  from  the  adsorber  during  the  first  3 
cycles  in  Figure  51,  the  maximal  temperature  in  the 
rate  in  Figure  53.  These  figures  show  that  the  qu 
second  cycle. 


Figure  51.  The  outlet  eoneentration  during  TSA  in 
the  1-eartridge  adsorber  (3  eyeles) 


adsorbent  bed  in  Figure  52  and  the  condensation 
asi-steady  state  is  practically  achieved  after  the 


Figure  52.  The  maximal  adsorbent  temperature 
during  TSA  in  the  1-eartridge  adsorber  (3  eyeles) 


Figure  53.  The  eondensation  rate  during  TSA  in  the  1- 
eartridge  adsorber  (3  eyeles) 
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In  Figure  53  an  enlarged  segment,  showing  the  desorption  step  during  the  seeond  eyele,  is  also  given. 
From  this  segment,  it  ean  be  seen  that  at  the  end  of  the  desorption  step,  the  eondensation  hasn’t  been 
finished. 


4.2.  Simulation  results  for  the  2-cartridges  adsorber 
4.2,1.  Adsorption  in  the  2-cartridges  adsorber 

Owing  to  larger  amount  of  adsorbent,  the  time  needed  to  eomplete  the  adsorption  proeess  in  the  2- 
eartridges  adsorber  is  mueh  longer  than  for  the  1-eartridge  adsorber.  For  that  reason,  these  simulations 
were  run  for  11000  s  or  ~  3  hours.  The  results  presented  here  were  obtained  for  the  same  inlet 
eoneentration  and  flow-rate  as  for  the  1-eartridge  adsorber  (0.001  mol/mol  and  0.1  mol/s). 

Velocity  distribution 

The  veloeity  field  in  the  2-eartridges  adsorber  is  presented  in  Figures  54-57.  In  Figure  54,  the  veloeity 
distribution  at  the  end  of  adsorption  (at  110000  s)  is  shown  in  the  form  of  a  2-D  surfaee  plot, 
together  with  a  veetor  representation  of  the  veloeity.  The  veloeity  distribution  is  also  presented  in  1-D 
diagrams:  the  time  profiles  at  the  middle  eross-seetions  of  both  eartridges  (z=16  em  and  z=42  em)  and 
4  radial  positions  in  Figure  55,  the  radial  profiles  at  z=16  em  and  z=42  em  in  Figure  56  and  the  axial 
profiles  in  several  radial  positions  in  Figure  57. 
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Figure  54.  Velocity  distribution  of  the 
gas  phase  during  adsorption  in  the  2- 
cartridge  adsorber  for  t=\  1000s 


Figure  55.  Time  profiles  of  the  gas  velocity  during 
adsorption  in  the  2-cartridge  adsorber  at  z=  16  cm 
(solid  lines)  and  z=49  cm  (dashed  lines) 
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Figure  56.  Radial  profiles  of  the  gas  veloeity 
during  adsorption  in  the  2-eartridge  adsorber  for 
^11000  s 


Figure  57.  Axial  profiles  of  the  gas  veloeity  during 
adsorption  in  the  2-eartridge  adsorber  at  r=1.25  em 


Concentration  distributions 

The  2-D  surface  plots  of  the  gas  and  solid  concentration  distributions  during  adsorption  in  the  2- 
cartridges  adsorber  are  given  in  Figures  58  and  59,  respectively.  The  concentrations  at  three  distinct 
moments  (^=100,  400  and  11000  s)  are  shown.  For  11000  s  the  gas  concentration  in  the  adsorber  is 
uniform  and  equal  to  the  feed  concentration  (0.001  mol/mol),  i.e.,  the  adsorption  process  is  finished. 
At  early  stages  of  adsorption,  the  concentrations  in  the  adsorber  are  nonuniform. 

The  concentration  profiles  are  also  given  in  1-D  diagrams:  the  time  profiles  at  z=  16  cm  in  Figure  60, 
the  radial  profiles  at  z=16  cm  in  Figure  61  and  the  axial  profiles  for  ^=200  s  in  Figure  62.  Figures  60a, 
61a  and  62a  correspond  to  the  gas  phase  and  Figures  60b,  61b  and  62b  to  the  solid  phase. 
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Figure  60.  Gas  concentration  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  /=100s,  (b)  /=400s  and 
(c)/=l  1000s 
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Figure  61.  Solid  concentration  distribution  during  adsorption  in  the  2-cartridges  adsorber:  (a)  /=100s,  (b)  /=400s  and 
(c)/=l  1000s 
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Figure  62.  Time  profiles  of  the  concentrations  in  the  2-cartridges  adsorber,  at  z=16cm  (solid  lines)  and  at  z=49cm 
(dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  63.  Radial  profiles  of  the  concentrations  during  adsorption  in  the  2-cartridges  adsorber,  at  z=  16cm  (solid 
lines)  and  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  64.  Axial  profiles  of  the  eoneentrations  during  adsorption  in  the  2-eartridges  adsorber,  for  /=200  s: 
(a)  gas  phase;  (b)  solid  phase 


Considerable  differenees  between  the  eoneentrations  in  the  two  eartridges  ean  be  observed.  These 
differenees  diminish  in  time,  as  both  eartridges  beeome  nearly  saturated. 


Temperature  distributions 

The  gas  and  solid  temperature  distributions  during  adsorption  in  the  2-eartridges  adsorber  are  given  in 
Figures  65  and  66,  respeetively,  in  the  form  of  2-D  surfaee  plots.  The  simulated  temperatures  are  also 
shown  in  1-D  diagrams:  the  time  profiles  in  Figures  67a  and  67b,  the  radial  profiles  in  Figures  68a  and 
68b  and  the  axial  profiles  in  Figures  69a  and  69b.  All  these  figures  show  temperature 
nonhomogeneity  in  both  phases.  The  maximal  temperature  rise  during  adsorption  is  less  than  4  ®C,  for 
the  2-eartridges  adsorber.  The  highest  temperature  is  observed  after  approximately  400  s  from  the  start 
of  adsorption. 
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Figure  65.  Gas  temperature  distribution  during  adsorption  in  the  2-eartridges  adsorber:  (a)  /=100s,  (b)  ^=400s  and 
(e)/=l  1000s 
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Figure  66.  Solid  temperature  distribution  during  adsorption  in  the  2-eartridges  adsorber:  (a)  /=100s,  (b)  /=400s  and 
(e)/=11000s 

(a)  (b) 


Figure  67.  Time  profiles  of  the  temperatures  during  adsorption  in  the  2-eartridges  adsorber,  atz=16em  (solid  lines) 
and  z=49em  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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r  [cm]  r  [cm] 

Figure  68.  Radial  profiles  of  the  temperatures  during  adsorption  in  the  2-eartridges  adsorber,  atz=16em  (solid 
lines)  and  z=49em  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  69.  Axial  profiles  of  the  temperatures  during  adsorption  in  the  2-eartridges  adsorber,  for  /=200  s: 
(a)  gas  phase;  (b)  solid  phase 

Large  differences  between  the  temperature  profiles  in  the  two  cartridges  can  also  be  observed. 


4.2.2.  Desorption  in  the  2-cartridges  adsorber 

Again,  the  complete  desorption  step  is  simulated,  by  using  a  combination  of  Model_D_2  and 
Model_DC_2.  The  results  presented  here  were  obtained  for  initial  concentration  in  the  solid  phase 
0.00465  mol/g,  electric  voltage  11  V  and  inert  gas  flow-rate  0.01  mol/s  per  cartridge  (0.002  mol/s  in 
total).  The  time  scale  used  for  this  simulation  is  the  same  as  for  the  1 -cartridge  adsorber:  0  to  200s. 

Velocity  distribution 

The  velocity  distribution  at  the  end  of  desorption  (at  ^=200  s)  is  shown  in  Figure  70,  in  the  form  of  a  2- 
D  surface  plot,  together  with  vector  representation  of  the  velocity.  The  velocity  distribution  is  also 
presented  in  1-D  diagrams:  the  time  profiles  in  Figure  71,  the  radial  profiles  in  Figure  72  and  the  axial 
profiles  in  Figure  73.  Some  (not  significant)  change  of  the  velocity  profile  in  time  is  present,  as  a 
result  of  the  change  of  the  gas  temperature. 


Figure  70.  Velocity  distribution  of  the 
gas  phase  during  desorption  in  the  2- 
cartridge  adsorber  at  /=200  s 


Figure  71.  Time  profiles  of  the  gas  velocity  during 
desorption  in  the  2-cartridge  adsorber,  at  z=  16cm 
(solid  lines)  and  at  z=49cm  (dashed  lines) 
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Figure  72.  Radial  profiles  of  the  gas  veloeity 
during  desorption  in  the  2-eartridge  adsorber,  for 
/=200  s 


Figure  73.  Axial  profiles  of  the  gas  veloeity  during 
desorption  in  the  2-eartridge  adsorber,  for  t=  200  s 


Concentration  distributions 

The  gas  and  solid  concentration  distributions  during  desorption  in  the  2-cartridges  adsorber,  presented 
in  the  form  of  2-D  surface  plots,  are  given  in  Figures  74  and  75,  respectively.  The  concentrations  at 
three  distinct  moments  (^=30,  90  and  200  s)  are  shown  (^30  s  is  the  moment  when  condensation 
starts,  i.e.  when  the  program  switches  from  Model_D_2  to  Model_DC_2).  The  concentration  profiles 
are  also  shown  in  1-D  diagrams:  the  time  profiles  in  Figure  76,  the  radial  profiles  in  Figure  77  and  the 
axial  profiles  in  Figure  78.  Figures  76a  to  78a  correspond  to  the  gas  phase  and  Figures  76b  to  77b  to 
the  solid  phase. 
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Figure  74.  Gas  concentration  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and  (c) 
/=200s 


Contrary  to  the  case  of  adsorption,  for  desorption  the  concentration  distributions  in  the  two  cartridges 
are  similar.  This  is  a  result  of  the  symmetrical  inlet  streams  in  the  two  cartridges. 
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Figure  75.  Solid  concentration  distribution  during  desorption  in  the  2-cartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and  (c)  /=200s 
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Figure  76.  Time  profiles  of  the  concentrations  during  desorption  in  the  2-cartridge  adsorber,  atz=16cm  (solid 
lines)  and  at  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  77.  Radial  profiles  of  the  concentrations  during  desorption  in  the  2-cartridge  adsorber,  at  z=  16cm  (solid 
lines)  and  at  z=49cm  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 
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Figure  78.  Axial  profiles  of  the  eoneentrations  during  desorption  in  the  2-eartridge  adsorber,  for  /=80  s: 
(a)  gas  phase;  (b)  solid  phase 


Condensation  flux  and  condensation  rate 

The  condensation  flux,  calculated  by  using  equation  (52)  is  shown  in  Figure  79  (the  time  profiles)  and 
Figure  80  (the  axial  profiles). 


Figure  79.  Time  profiles  of  the  eondensation  flux  in  Figure  80.  Axial  profiles  of  the  eondensation  flux  in 

the  2-eartridge  adsorber  2-eartridges  adsorber 
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Figure  81.  Condensation  rate  in  the  2-eartridges  adsorber 


These  figures  show  that  condensation  starts  at 
^=30  s  and  the  condensation  flux  is  maximal  at 
approximately  f=60s.  At  ^=200s,  when  the 
power  is  switched  of,  the  condensation  is  still 
not  completely  finished.  At  the  beginning  of  the 
process  the  condensation  flux  increases  from 
the  adsorber  ends  towards  its  middle,  with  a 
sharp  minimum  at  the  middle  cross-section, 
between  the  two-cartridges. 

Figure  81  shows  the  time  profile  of  the 
condensation  rate,  calculated  by  integration  of 
the  condensation  flux  along  the  adsorber  wall 
surface. 
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Temperature  distributions 

The  gas  and  solid  temperature  distributions  in  the  2-eartridges  adsorber  in  whieh  desorption  is  taking 
plaee  are  again  presented  in  the  form  of  2-D  surfaee  plots  (Figures  82  and  83)  and  in  1-D  plots:  the 
time  profiles  in  Figures  84a  and  84b,  the  radial  profiles  in  Figures  85a  and  85b  and  the  axial  profiles  in 
Figures  86a  and  86b.  All  these  figures  show  temperature  nonhomogeneity  in  both  phases  and  large 
temperature  ehanges  during  desorption.  Similarity  between  the  temperature  distributions  in  the  two 
eartridges  ean  easily  be  observed. 
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Figure  82.  Gas  temperature  distribution  during  desorption  in  the  2-eartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and  (e)  /=200s 
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Figure  83.  Solid  temperature  distribution  during  desorption  in  the  2-eartridges  adsorber:  (a)  /=30s,  (b)  /=90s  and  (e)  /=200s 
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Figure  84.  Time  profiles  of  the  temperatures  during  desorption  in  the  2-eartridges  adsorber,  atz=16em  (solid  lines) 
and  z=49em  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  85.  Radial  profiles  of  the  temperatures  during  desorption  in  the  2-eartridges  adsorber,  at  z=16em  (solid 
lines)  and  z=49em  (dashed  lines):  (a)  gas  phase;  (b)  solid  phase 


(a)  (b) 


Figure  86.  Axial  profiles  of  the  temperatures  during  desorption  in  the  2-eartridges  adsorber,  for  /=80  s 
(a)  gas  phase;  (b)  solid  phase 
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Electric  resistivity  [Qcm] 


Distribution  of  the  electric  resistivity  of  the  ACFC  adsorbent  bed 

The  distribution  of  the  eleetrie  resistivity  in  the  2-eartridges  adsorber  is  shown  in  Figure  87  (in  the 
form  of  2-D  surfaee  plot)  and  in  Figures  88,  89  and  90  (the  1-D  diagrams  showing  the  time,  radial  and 
axial  distributions,  respeetively). 
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Figure  87.  Distribution  of  the  eleetrie  resistivity  of  the  ACFC  adsorbent  in  the  2-eartridges  adsorber:  (a)  /=30s,  (b)  ^90s  and 
(e)  /=200s 


Figure  88.  Time  profiles  of  the  eleetrie  resistivity  in 
the  2-eartridges  adsorber,  atz=16em  (solid  lines)  and 
at  z=49em  f dashed  lines') 


Figure  89.  Radial  profiles  of  the  eleetrie  resistivity  in 
the  2-eartridges  adsorber,  atz=16em  (solid  lines)  and 
z=49em  f dashed  lines') 


The  time  dependenee  of  the  eleetrie  power  is  shown  in  Figure  91. 
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Maximal  adsorbent  temperature  [K]  Outlet  concentrtion  [moi/moi] 


Figure  90.  Axial  profiles  of  the  eleetrie  resistivity  in  Figure  91.  Used  eleetrie  energy  during  desorption  in  the  2- 

the  2-eartridges  adsorber,  for  ^80  s  eartridges  adsorber 


4.2.3.  TSA  process  in  the  2-cartridges  adsorber 


Figure  92.  Time  profiles  of  the  outlet  eoneentration 
in  the  2-eartridges  adsorber,  for  3  eyeles 


The  simulation  results  presented  here  were 
obtained  for  7^^=463.15  K  and  Cbreak=0.5  Qn. 
All  other  simulation  parameters  are  the  same  as 
in  the  simulations  of  adsorption  and  desorption. 
The  ehange  of  the  outlet  eoneentration  from  the 
adsorber  during  the  first  3  eyeles  in  is  shown  in 
Figure  92,  the  maximal  temperature  in  the 
adsorbent  beds  in  Figure  93  and  the 
eondensation  rate  in  Figure  94.  The  quasi¬ 
steady  state  is  praetieally  aehieved  after  the 
seeond  eyele. 

It  should  be  notieed  that  the  TSA  eyele  in  the  2- 
eartridges  adsorber  lasts  eonsiderably  longer 
that  in  the  1-eartridge  adsorber.  This  is 
primarily  a  result  of  a  mueh  longer  adsorption 
half-eyele. 


Figure  93.  Maximal  adsorbent  temperature  in  the 
2-eartridges  adsorber  for  3  eyeles 
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Figure  94.  Condensation  rate  in  the  2-eartridges 
adsorber  for  3  eyeles 
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5.  Analysis  of  the  influence  of  process  parameters  and  optimization  of  the  TSA  process 

In  order  to  analyze  the  results  of  the  TSA  proeesses,  some  eriterial  parameters  should  be  defined,  as 
measures  of  the  proeess  quality.  We  define  the  following  eriterial  parameters: 

•  Adsorption  time  (duration  of  the  adsorption  half-eyele),  ta 

•  Desorption  time  (duration  of  the  desorption  half-eyele),  Td 

•  Average  outlet  eoneentration  during  adsorption, 

•  Average  outlet  eoneentration  during  desorption,  <Com^>d 

•  Separation  faetor  defined  as  the  ratio  of  the  average  outlet  eoneentrations  during  adsorption 
and  desorption,  ^s-^Cour^Ti!  <Cout>A 

•  Purifieation  faetor,  defined  as  the  ratio  of  the  outlet  and  inlet  eoneentrations  during  adsorption, 

a!  ^in,A 

•  The  total  amount  of  eondensed  liquid  per  eyele,  Lcond 

•  Used  eleetrie  energy  per  eyele,  Qei 

•  Used  energy  per  mol  of  eondensate,  QeilLcond 

In  prineiple,  the  following  ean  be  eonsidered  as  favourable: 

Long  adsorption  time 

Short  desorption  times 

Low  outlet  eoneentration  during  adsorption 

High  outlet  eoneentrations  during  desorption 

High  separation  faetor 

Low  purifieation  faetor 

Large  amount  of  eondensed  liquid  per  eyele 

Low  energy  spent  per  eyele 

Low  energy  spent  per  mole  of  eondensate. 

In  this  seetion,  we  will  analyze  how  these  eriteria  are  influeneed  by  the  ehange  of  some  input 
parameters  of  the  TSA  proeess.  The  influenee  of  the  following  parameters  was  ehosen  for  analysis: 

The  maximal  adsorbent  temperature  at  whieh  the  TSA  eyele  is  switehed  from  desorption 
to  adsorption  -  the  switeh  temperature 

The  breakthrough  eoneentration  for  whieh  the  TSA  eyele  is  switehed  from  adsorption  to 
desorption,  Cbreah 

The  eleetrieal  voltage  used  during  desorption  Uq; 

The  inert  flow-rate  during  the  adsorption  part  of  the  eyele  Ga; 

The  inert  flow-rate  during  the  desorption  part  of  the  eyele  Gd. 

The  analysis  was  performed  by  ehanging  one  of  these  five  parameters  at  a  time  and  keeping  all  the 
others  eonstant.  The  “fixed”  parameters  have  the  same  values  as  in  the  simulations  presented  in  the 
previous  seetion.  For  eaeh  eombination  of  parameters  the  9  eriteria  defined  previously  were 
ealeulated.  The  results  are  presented  in  tabular  form. 

The  results  presented  in  this  seetion  eorrespond  to  the  seeond  eyeles  of  the  simulated  TSA  proeesses. 
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The  results  obtained  for  the  1-eartridge  adsorber  are  presented  first. 


5.1.  Analysis  for  the  1-cartridge  adsorber 

In  Table  3  we  present  the  eriterial  parameters  ehosen  for  analysis  obtained  for  three  different  switeh 
temperatures  (the  maximal  temperature  of  the  adsorbent  bed  at  whieh  the  TSA  proeess  is  switehed 
from  adsorption  to  desorption). 


Table  3.  Analysis  of  the  influenee  of  the  switeh  temperature,  for  the  1-eartridge  adsorber 


Criteria 

Switch  temperature  Tsw  [K] 

433.15 

463.15 

493.15 

<C«f>D  [mol/mol] 

0.00853 

0.00879 

0.00895 

<Coui>K  [mol/mol] 

2.65e-04 

2.00e-04 

1 .84e-04 

32.189 

43.950 

48.641 

A  /^in,A 

0.265 

0.200 

0.184 

Lcond  [mol] 

0.0400 

0.0499 

0.0557 

Qel  [KJ] 

13.363 

16.735 

19.593 

Qel  !  ^cond  [KJ/lHOl] 

334.07 

335.37 

351.76 

Td  [s] 

130 

158 

180 

Ta  [s] 

2370 

2440 

2600 

These  results  show  that  the  following  eriteria  inerease  with  inereasing  the  switeh  temperature:  the 
adsorption  and  desorption  times,  the  outlet  eoneentration  during  desorption  and  the  separation  faetor, 
the  eondensed  liquid  and  the  used  energy  per  eyele.  On  the  other  hand,  the  outlet  eoneentration  during 
adsorption  and  the  purifieation  faetor  deerease  with  inereasing  The  energy  eonsumption  per  mole 
of  eondensate  is  nearly  eonstant.  It  eould  be  eoneluded  that  the  inerease  of  the  switeh  temperature  is 
favorable  from  the  point  of  view  of  separation  and  purifieation  of  the  feed  stream,  as  well  as 
regeneration  of  the  adsorbed  vapor,  but  unfavorable  from  the  point  of  view  of  energy  eonsumption.  It 
should  also  be  remembered,  that  the  switeh  temperature  is  often  limited  by  the  temperature  of 
degradation  of  some  thermally  unstable  materials  used  in  the  proeess. 

In  Table  4  we  show  the  influenee  of  the  breakthrough  eoneentration,  i.e.,  the  outlet  gas  eoneentration 
for  whieh  the  TSA  proeess  is  switehed  from  adsorption  to  sorption. 


Table  4.  Analysis  of  the  influenee  of  the  breakthrough  eoneentration  Ctreah  for  the  1-eartridge  adsorber 


Criteria 

Breakthrouj 

concentration  Cfyreak  [mol/mol]  | 

2.5e-4  (25%  C,„) 

3.5e-4  (35%  C,„) 

5.0e-4  (50%  C,„) 

<Coui>D  [mol/mol] 

0.00835 

0.00854 

0.00879 

<Coui>A  [mol/mol] 

1 .28e-04 

1 .58e-04 

2.00e-04 

65.234 

54.051 

43.950 

A  /Cin,A 

0.128 

0.158 

0.200 

^cond  [mol] 

0.0397 

0.0437 

0.0499 

Qel  [KJ] 

15.268 

15.717 

16.735 

Qel  !  ^cond  [KJ/lHOl] 

384.58 

359.Q6 

335.37 

Td  [s] 

143 

148 

158 

Ta  [s] 

1800 

2050 

2440 
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The  results  given  in  Table  4  show  that  the  inerease  of  the  breakthrough  eoneentration  results  in  longer 
adsorption  and  desorption  times,  higher  outlet  eoneentrations  during  both  adsorption  and  desorption, 
lower  separation  and  higher  purifieation  faetors,  larger  amounts  of  eondensed  liquid  and  energy  spent 
per  eyele,  but  lower  energy  eonsumption  per  mole  of  eondensed  liquid. 

Table  5  shows  the  eriterial  parameters  eorresponding  to  simulations  obtained  with  3  different  values  of 
the  supply  voltage  during  desorption. 


Table  5.  Analysis  of  the  influenee  of  the  voltage  used  for  desorption,  f/o,  for  the  1-eartridge  adsorber 


Criteria 

Supply  voltage  Uq  [V]  | 

10 

11 

12 

<C«f>D  [mol/mol] 

0.00895 

0.00879 

0.00869 

<C„ui>K  [mol/mol] 

1 .96e-04 

2.00e-04 

2.08e-04 

45.663 

43.950 

41.779 

A  /^in,A 

0.196 

0.200 

0.208 

Lcond  [mol] 

0.0500 

0.0499 

0.0510 

Qe,  [KJ] 

19.626 

16.735 

15.433 

Qel  /  Lcond  [KJ/mol] 

392.52 

335.37 

302.61 

Td  |s] 

224 

158 

122 

Ta  [s] 

2460 

2440 

2430 

When  the  supply  voltage  during  desorption  inereases  the  slightly  shorter  adsorption  and  eonsiderably 
shorter  desorption  times  are  obtained.  The  outlet  eoneentrations  during  adsorption  and  desorption  both 
inerease  with  the  inerease  of  f/o,  but  only  slightly,  as  well  as  the  purifieation  faetor.  The  separation 
faetor  deereases  somewhat,  and  the  used  energy  per  eyele  and  per  mole  of  eondensate  eonsiderably. 
The  amount  of  eondensed  liquid  per  eyele  doesn’t  ehange  eonsiderably. 

Tables  6  shows  the  influenee  of  the  flow-rate  during  adsorption,  GA,  and  Table  7  of  the  flow-rate 
during  desorption,  GD,  on  the  eriteria  ehosen  for  eomparison.  For  the  ehosen  range  of  flow-rates  and 
for  the  parameters  used  for  simulation,  these  influenees  are  not  signifieant. 


Table  6.  Analysis  of  the  influenee  of  the  flow-rate  during  adsorption  Ga,  for  the  1-eartridge  adsorber 


Criteria 

Flow-rate  during  adsorption  Ga  [mol/s]  | 

0.08 

0.09 

0.10 

<Coui>u  [mol/mol] 

0.00878 

0.0088 

0.00879 

<Coui>A  [mol/mol] 

2.02E-04 

2.01  E-04 

2.00e-04 

OCs  ^^ouP^D  l^^ouP^ A 

44.7959 

43.7811 

43.950 

^p~^^ouP^ A  /Cin,A 

0.202 

0.201 

0.200 

Lcond  [niol] 

0.0502 

0.0506 

0.0499 

Qel  [KJ] 

16.737 

16.842 

16.735 

Qel !  LcQfifi  [KJ/mol] 

333.41 

332.85 

335.37 

Td  [s] 

158 

159 

158 

Ta  [s] 

2470 

2460 

2440 
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Table  7.  Analysis  of  the  influenee  of  the  flow-rate  during  desorption  Gd,  for  the  1-eartridge  adsorber 


Criteria 

Flow-rate  during  desorption  Gd  [mol/s]  | 

0.01 

0.02 

0.03 

<C„ui>o  [mol/mol] 

0.00879 

0.00868 

0.00886 

<C„ui>K  [mol/mol] 

2.00e-04 

2.00e-04 

2.00e-04 

43.950 

44.8731 

43.9801 

0.200 

0.200 

0.200 

Lcond  [mol] 

0.0499 

0.0528 

0.0526 

Qel  [KJ] 

16.735 

17.471 

17.900 

Qel  /  Lcond  [KJ/mol] 

335.37 

330.89 

340.30 

Td  [s] 

158 

165 

169 

Ta  [s] 

2440 

2410 

2430 

5.2.  Analysis  for  the  2-cartridges  adsorber 

As  our  analysis  for  the  1-eartridge  adsorber  showed  that  the  flow-rates  don’t  influenee  the  TSA 
proeess  mueh,  only  three  input  parameters  were  varied  for  the  2-eartridges  adsorber:  the  switeh 
temperature,  the  breakthrough  eoneentration  and  the  eleetrie  voltage.  Their  influenee  on  the  analyzed 
eriterial  parameters  of  he  TSA  proeess  is  shown  in  Tables  8,  9  and  10,  respeetively. 


Table  8.  Analysis  of  the  influenee  of  the  switeh  temperature  Tsw,  for  the  2-eartridges  adsorber 


Criteria 

Switch  temperature  Tsw  [K]  | 

433.15 

463.15 

485.28 

<C„ui>M  [mol/mol] 

0.00921 

0.00933 

0.00936 

<C„ui>K  [mol/mol] 

2.78E-04 

2.30E-04 

2.22E-04 

Ots  l^^ouP^ A 

33.1295 

40.4948 

42.1622 

^p~^^ouP^ A  /Cin,A 

0.278 

0.230 

0.222 

Lcond  [niol] 

0.114 

0.1270 

0.131 

Qel  [KJ] 

34.648 

41.890 

44.665 

Qel !  LcQfifi  [KJ/mol] 

303.93 

329.84 

340.95 

Td  [s] 

168 

198 

208 

Ta  [s] 

8090 

8570 

8720 

Table  9.  Analysis  of  the  influenee  of  the  breakthrough  eoneentration  C break,  for  the  2-eartridges  adsorber 


Criteria 

Breakthrouj 

»h  concentration  Ctreak  [mol/mol] 

2.5e-4  (25%  C,„) 

3.5e-4  (35%  C,„) 

5.0e-4  (50%  C,„) 

<C«f>D  [mol/mol] 

0.00914 

0.00917 

0.00933 

<Coui>K  [mol/mol] 

1 .49E-04 

1 .92E-04 

2.30E-04 

ouP^H  /"^L^ouP^ A 

61 .3423 

47.7604 

40.4948 

^L^ouP^ A  /^in,A 

0.149 

0.192 

0.230 

Lcond  [mol] 

0.0879 

0.1030 

0.1270 

Qel  [KJ] 

35.986 

37.531 

41 .890 

Qel  /  Lcond  [KJ/mol] 

409.40 

364.38 

329.84 

Td  [s] 

170 

179 

198 

Ta  [s] 

5720 

6840 

8570 
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Table  10.  Analysis  of  the  influenee  of  the  voltage  supply  f/o,  for  the  2-eartridges  adsorber 


Criteria 

Voltage  supply  t/o  [V]  | 

10 

11 

12 

<C„ui>M  [mol/mol] 

0.00937 

0.00933 

0.00924 

<Ci,f>A  [mol/mol] 

2.29E-04 

2.30E-04 

2.51  E-04 

40.9170 

40.4948 

36.8127 

0.229 

0.230 

0.251 

Lcond  [mol] 

0.1197 

0.1270 

0.1235 

Qel  [KJ] 

45.837 

41.890 

35.899 

Qel  !  ^cond  [KJ/mol] 

382.93 

329.84 

290.68 

Td  [s] 

265 

198 

143 

Ta  [s] 

8440 

8570 

8600 

Similar  eonelusions  as  for  the  1-eartridge  adsorber  ean  be  withdrawn. 

Comparison  of  the  results  obtained  for  the  1-eartridge  and  2-eartridges  adsorbers  shows  that  the 
adsorption  times  are  more  than  double  and  the  desorption  times  about  50%  longer  for  the  2-eartridges 
adsorber.  The  separation  and  the  purifieation  faetors  are  similar  for  both  adsorbers.  On  the  other  hand, 
the  amount  of  eondensed  liquid  and  the  used  energy  per  eyele  for  the  2-eartridges  adsorber  are  more 
than  double  than  for  the  1-eartridge  adsorber.  Nevertheless,  the  values  of  the  energy  spent  per  mole  of 
eondensate  are  similar  for  both  adsorber  types. 


5.3.  Discussion  of  the  results  and  some  recommendations  for  choosing  the  optimal  input 
parameters 

The  inspeetion  of  the  results  listed  in  Tables  3-10  shows  that  unexpeetedly  low  values  of  the  amount 
of  eondensed  liquid  per  eyele  are  obtained.  We  found  the  reason  for  that  in  the  faet  that  in  our  model 
the  desorption  step  is  ended  when  the  predefined  switeh  temperature  is  reaehed,  while  the 
eondensation  proeess  is  still  not  finished,  as  ean  be  seen  from  Figure  53.  After  that  moment,  the 
simulation  switehes  to  the  model  for  adsorption  whieh  doesn’t  assume  eondensation.  In  order  to 
resolve  this  impreeision,  an  additional  segment  should  be  ineorporated  in  the  integral  TSA  model,  for 
an  additional  step  of  the  TSA  proeess  of  desorption  with  in- vessel  eondensation,  but  without  eleetrieal 
heating.  We  expeet  that  sueh  a  model  would  prediet  larger,  more  realistie  amounts  of  eondensate  per 
eyele,  espeeially  for  the  eases  with  short  desorption  times.  The  spent  energy  per  eyele  would  remain 
unehanged,  so  we  expeet  that  mueh  lower  values  of  the  energy  eonsumption  per  mole  of  eondensate 
would  be  obtained. 

Regarding  the  influenee  of  the  input  parameters  on  the  TSA  eriterial  parameters,  there  is  no  ideal 
solution,  i.e.  all  input  parameters  influenee  some  of  the  eriteria  in  a  favorable  and  some  in  an 
unfavorable  way.  For  that  reason,  the  reeommendations  for  ehoosing  optimal  input  parameters  depend 
on  the  priorities  of  the  eriterial  parameters.  Here  are  some  examples: 

If  purification  of  the  gas  stream  is  the  main  priority,  lower  values  of  the  breakthrough 
eoneentration  should  be  used,  primarely.  Higher  switeh  temperature  and  lower  supply 
voltage,  also  result  with  more  effieient  purifieation,  to  some  extent. 

If  regeneration  of  the  vapor  in  the  liquid  form  is  the  priority  of  the  proeess,  one  should 
use  high  switeh  temperature,  large  breakthrough  eoneentration  and  high  voltage. 

If  energy  efficiency  of  the  proeess  is  the  primary  goal,  low  switeh  temperature,  large 
breakthrough  eoneentration  and  high  voltage  should  be  used. 
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6.  Conclusions  and  possibilities  for  model  improvement 


In  this  phase  of  the  projeet  we  have  developed  two  eomplex  models  for  simulation  of  integral  TSA 
eyeles  in  adsorbers  with  one  or  two  annular,  radial-flow,  eartridge-type  adsorption  beds.  Eaeh  model 
was  built  of  three  segments,  eorresponding  to  three  steps  of  the  TSA  proeess:  adsorption, 
eleetrothermal  desorption  and  eleetrothermal  desorption  with  in-vessel  eondensation.  The  individual 
segments  were  built  using  Femlab  software,  while  Femlab  with  Matlab  was  used  for  building  the 
integral  TSA  models. 

The  models  were  sueeessfully  used  for  simulation  of  separate  stages  and  eomplete  TSA  eyeles  in  the 
1-eartridge  and  2-eartridges  adsorbers  and  for  analysis  of  the  influenee  of  the  main  input  parameters 
on  the  system  performanee.  Some  basie  reeommendations  regarding  the  ehoiee  of  the  input  parameters 
of  the  TSA  proeess  were  given  based  on  it,  in  order  to  show  that  the  developed  models  ean  be  used  for 
optimization. 

The  models  were  developed  and  run  on  Pentium  4  PCs,  with  1  GB  RAM  memory.  During  their 
development,  eonvergenee  problems  oeeurred  frequently.  Regarding  eonvergenee,  the  models  are 
pretty  sensitive  to  the  ehoiee  of  the  model  parameters.  The  eomputer  times  needed  for  simulations 
were  generally  long,  but  differed  from  ease  to  ease.  In  prineiple,  the  eomputer  times  needed  to  solve 
the  2-eartridges  model  were  mueh  longer  than  for  the  1-eartridge  model. 


In  spite  of  the  great  eomplexity  of  the  developed  models,  there  are  still  some  possibilities  for  their 
upgrading  in  order  to  make  them  more  realistie.  Some  of  them  would  be: 

Adding  a  model  segment  defining  the  end  of  the  desorption  half-eyele  without  heating, 
between  switehing  the  power  off  and  the  end  of  eondensation. 

Adding  relations  for  ealeulation  of  some  physieal  and  transport  parameters  and 
eoeffieients  as  fiinetions  of  the  temperatures,  eoneentrations  and  veloeities  in  the  system, 
instead  of  using  eonstant  values. 

Adding  a  model  for  the  heat  balanee  of  the  adsorber  wall. 
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Nomenclature 


a 

C  (mol/mol) 
C*(mol/mol) 
Cbreak  (mol/mol) 
Csat  (mol/mol) 
Cpg  (J/mol/K) 

Cpi  (J/mol/K) 

Cps  (J/g/K) 

Cpv  (J/mol/K) 
ZJnir  (mol/cm/s) 
Dmz  (mol/cm/s) 
a’'®  (W/K/cm) 
A***  (W/K/cm) 

E  (J/mol) 
g  (cm/s^) 

G  (mol/s) 

H  {cm) 
hi  (cm) 
hi  (cm) 
h  (J/cm^/K) 

Ki  (J/cm^/K) 

/Zs2  (J/cm^/K) 

Awg  (J/cm^/K) 

J  (A/cm^) 

Jcond  (mol/cm^/s) 
k  (cm^) 

km  (mol/cm^/s) 

km\  (mol/cm^/s) 

kmi  (mol/cm^/s) 
tube 


-  Specific  surface  area 

-  Temperature  coefficient  of  the  bed  electrical  resistivity 

-  Adsorbate  concentration  in  the  gas  phase 

-  Adsorbate  concentration  in  the  gas  phase  in  equilibrium  with  the  solid  phase 

-  Breakthrough  concentration 

-  Saturation  concentration 

-  Specific  heat  capacity  of  the  inert  gas 

-  Specific  heat  capacity  of  liquid  adsorbate 

-  Heat  capacity  of  the  solid  phase 

-  Specific  heat  capacity  of  gaseous  adsorbate 

-  Radial  mass  transfer  dispersion  coefficient 

-  Axial  mass  transfer  dispersion  coefficient 

-  Heat  diffiisivity  of  the  gas  phase 

-  Heat  diffusivity  of  the  solid  phase 

-  Adsorption  energy  of  the  adsorbate  (D-R  eq.) 

-  Gravitation  constant 

-  Flow  rate  of  the  inert  gas 

-  Bed  axial  dimension  (Fig.  5) 

-  Adsorber  dimension  (Fig.  5) 

-  Adsorber  dimension  (Fig.  5) 

-  Solid  to  gas  heat  transfer  coefficient 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the 

-  Heat  transfer  coefficient  from  the  solid  phase  to  the 

-  Gas  to  ambient  heat  transfer  coefficient  (heat  losses) 

-  Current  density 

-  Condensation  flux 

-  Asorbent  bed  permeability 

-  Mass  transfer  coefficient  in  the  adsorbent  bed 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  central  tube 

-  Mass  transfer  coefficient  between  the  adsorbent  bed  and  the  gas  in  the  annular 


gas  phase  in  the  central  tube(s) 
gas  phase  in  the  annular  tube 


C«./(mol/s) 

l-‘cond  (mol) 

p(Pa) 

Pa  (Pa) 


-  Flow-rate  of  the  condensed  liquid 

-  Total  amount  of  the  condensed  liquid 

-  Gas  pressure 

-  Adsorbate  partial  pressure 
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Pc  (Pa) 

-  Critical  pressure 

/(Pa) 

-  Adsorbate  saturation  pressure 

a/(W) 

-  Rate  of  heat  generation  (equal  to  electric  power) 

QelH) 

-  Heat  generation  (equal  to  electric  energy) 

q  (mol/g) 

-  Adsorbate  concentration  in  the  solid  phase 

Rg  (J/mol/K) 

-  Gas  constant 

r  (cm) 

-  Radial  dimension 

ri  (cm) 

-  Adsorber  dimension  (Fig.  5) 

r2  (cm) 

-  Adsorber  dimension  (Fig.  5) 

rs  (cm) 

-  Adsorber  dimension  (Fig.  5) 

r4  (cm) 

-  Adsorber  dimension  (Fig.  5) 

7^a(K) 

-  Ambient  temperature 

7’g(K) 

-  Gas  phase  temperature 

7’c(K) 

-  Critical  temperature 

7’r(K) 

-  Referent  temperature 

n(K) 

-  Solid  phase  temperature 

7;w(K) 

-  Switch  temperature  (from  desorption  to  adsorption) 

-  Wall  temperature 

?(s) 

-  Time 

U(V) 

-  Electric  potential 

Uo(V) 

-  Supply  voltage 

u  (cm/s) 

-  Radial  superficial  gas  velocity 

VPa,  VPb,  VPc, 

VPd  -  Wagner  constants 

V  (cm/s) 

-  Axial  superficial  gas  velocity 

W  (cm^/g) 

-  Adsorbate  volume  per  1  g  of  the  adsorbent  mass 

Wo  (cmVg) 

-  Total  volume  of  micropores  (D-R  equation) 

Greek  letters: 

dp 

-  Purification  factor 

dr 

-  Regeneration  factor 

ds 

-  Separation  factor 

A/4ds  (J/mol) 

-  Heat  of  adsorption 

A//cond  (J/mol) 

-  Heat  of  condensation 

Eb 

-  Bed  porosity 

g  (Pa/s) 

-  Dynamic  viscosity 

p  (Qcm) 

-  Electric  resistivity 

Po  (Qcm) 

-  Electric  resistivity  at  referent  temperature  Tr 
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Pg  (mol/cm^) 

-  Gas  phase  density 

Pb  (g/cm^) 

-  Adsorbent  bed  density 

PA(g/cm^) 

-  Adsorbate  density 

x(s) 

-  ‘‘Half’  eyele 

Subscripts: 

A 

-  adsorption 

b 

-  bed 

D 

-  desorption 

g 

-  gas 

in 

-  inlet 

out 

-  outlet 

P 

-  previous  (initial) 

r 

-  radial 

s 

-  solid  phase 

it 

-  inner  (eentral)  tube 

ot 

-  outer  (annular)  tube 

z 

-  axial 

Special  characters: 

<>  -  time  average  value 
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